Colloquium: Nonlinear collective interactions in quantum plasmas with 
degenerate electron fluids 



P. K. Shukla* and B. Eliasson 

RUB International Chair, Fakultat fur Physik und Astronomie, Ruhr-Universitat Bochum, D-44780 Bochum, 
Germany 

(Received 27 September 2010; Revised 18 February 2011; Accepted 6 May 2011, in Reviews of Modern Physics) 

The current understanding of some important nonlinear collective processes in quantum plas- 
mas with degenerate electrons is presented. After reviewing the basic properties of quantum 
plasmas, we present model equations (e.g. the quantum hydrodynamic and effective nonlinear 
Schrodinger-Poisson equations) that describe collective nonlinear phenomena at nanoscales. The 
effects of the electron degeneracy arise due to Heisenberg's uncertainty principle and Pauli's exclu- 
sion principle for overlapping electron wavefunctions that result in tunneling of electrons and the 
electron degeneracy pressure. Since electrons are Fermions (spin-1/2 quantum particles), there 
also appears an electron spin current and a spin force acting on electrons due to the Bohr mag- 
netization. The quantum effects produce new aspects of electrostatic (ES) and electromagnetic 
(EM) waves in a quantum plasma that are summarized in here. Furthermore, we discuss non- 
linear features of ES ion waves and electron plasma oscillations (ESOs), as well as the trapping 
of intense EM waves in quantum electron density cavities. Specifically, simulation studies of the 
coupled nonlinear Schrodinger (NLS) and Poisson equations reveal the formation and dynamics 
of localized ES structures at nanoscales in a quantum plasma. We also discuss the effect of an 
external magnetic field on the plasma wave spectra and develop quantum magnctohydrodynamic 
(Q-MHD) equations. The results are useful for understanding numerous collective phenomena in 
quantum plasmas, such as those in compact astrophysical objects (e.g. the cores of white dwarf 
stars and giant planets), as well as in plasma-assisted nanotechnology (e.g. quantum diodes, 
quantum free-electron lasers, nanophotonics and nanoplasmonics, metallic nanostructures, thin 
metal films, semiconductor quantum wells and quantum dots, etc.), and in the next-generation 
of intense laser-solid density plasma interaction experiments relevant for fast ignition in incrtial 
confinement fusion schemes. 
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I. INTRODUCTION 

Dense plasmas composed of ions, degenerate electrons, 
positrons, and/or holes (in the context of semiconduc- 
tors) are referred to as quantum plasmas. In the latter, 
the degeneracy of the lighter plasma species (electrons, 
positrons, holes) appears at very high densities and rel- 
atively low temperatures, where the mean inter-particle 
distance is smaller than (or of the same order as) the de 
Broglie thermal wavelength. The ions are typically non- 
degenerate due to their relatively large mass in compari- 
son with the electron mass. In quantum physics, Heisen- 



[Tgl berg's uncertainty principle (Bransden and Joachain 



2000} | Dirac[ fl98Tj |Holland[ |1993[ |Landau and LifeTutz 
1998a|) dictates that conjugate variables, such as the po- 



sition and momentum of a particle cannot be precisely 
determined simultaneously; the product of the uncer- 



2 



tainties of the position and momentum is equal to or 
larger than h/2, where h (= 1.0544 x 10 -27 erg sec) is 
Planck's constant divided by 2ir. The position of an elec- 
tron subjected to the influence of an atomic nucleus is 
very well defined (the force to which it is subjected is 
large). However, owing to Heisenberg's uncertainty prin- 
ciple, the electron momentum is ill defined. An electron 
has a continuous motion around the position it occupies. 
This motion exerts pressure on the surrounding medium, 
exactly as the thermal agitation of the particles of a gas 
exerts its pressure. This pressure is called the electron de- 
generacy pressure. This pressure, since it is nonthermal 
in origin, is, of course, independent of the electron tem- 
perature; the pressure of degenerate electrons increases 
with increasing electron number density. It is, however, 
only at very high densities that the degeneracy pressure 
becomes comparable or larger to the thermal gas pres- 
sure. One then says that the plasma matter is in an ex- 
otic state, comprising degenerate electrons and positrons 
or holes. 

Plasmas with degenerate electrons and positrons with 
number densities comparable with solids and tempera- 
tures of several electron volts fall under the category 



powerful x-ray sources are employed for accessing narrow 
bandwidth spectral lines via collective Thomson scatter- 
ing of light off electron-density fluctuations. These exper- 
imental techniques also allow accurate measurements of 
the electron velocity distribution function, temperature, 
and ionization state in the dense matter regime. |Gre-| 
gori and Gericke ( 2009 1 also proposed future experiments 



to measure low-frequency oscillations in plasmas when 
keV free-electron lasers will become available. IFroula e~h 



al. (2011 1 summarizes the measurement techniques using 



scattering of EM waves in plasmas, and recent experi- 
mental results from x-ray scattering experiments in dense 
plasmas reveal that quantum mechanical effects are in- 



2007) 



dee d imp ortant (|Glenzer and Redmer 
al 



2009 Glenzer et 



Furthermore, due to recent experimental progress in 
femtosecond pump-probe spectroscopy, the field of quan- 



tum plasmas is also gaining significant attention (Crou 



seillcs et al. 20081 in connection with the collective dy- 



namics of an ensemble of degenerate electrons in metal- 
lic nanostructures and thin metal films. The physics of 
quantum plasmas is also relevant in the context of quan- 



of dense matter (Fortov 2009 Ichimaru 19821 that ap 



pears in the core of giant planets (Chabrier et al. 2006 



Chabrier 



(Guillot 



2009 
1999) 



son 



turn diodes (Ang et al. 2003 2007 Shukla and Elias- 



2008b), nanophotonics an d nanowires (|Barnes et 



HoTnl [19911 and the crusts of old stars nanop l asmonics QAtwater 



al.\ |2003| |Chang et al] |2006[ |Shpatakovskaya[ |2006|), 



Dense compressed plasmas are currently 



et al. 2008 Ozbay 2006f Stockman 2011), high-gain 



2007 Maier, 2007 Marklund 



of wide interest due to their applications to astrophysi- quantum free-electron l asers qSerbeto et a/.[|2008[|2009D , 



cal and cosmological environments (Benvenuto and De 
" [20051 IHarding and~Tai[ [20061 |Lai[ |2001t |Opher| 



Vito 



microplasma systems ( Becker et al. 2006 1 , and small 



2001), as well as to inertial fusion science involving in- 



Manfredi and Hervieux 



semiconductor devices (Haug and Koch 2004 2007 



2007 Markovich et al. 19901, 



tense laser-solid density plasma interaction experiments 



(| Andreev[ [20001 |Froula et al. \ |201 1 



Glenzer and Redmer 
2009[ jGlenzer et q/.[|2007[|Hu and keitel[ |1999||Kritcher| 
et al.\ 2008: Lee et al. 1 120091 ILindl. 1995: Malkin et al. 



such as quantum wells and piezomagnetic quantum dots 
(|Adolfatn et al.\ |2008h. The latter can be used 



as 



2007| IMarklund and Shukla] [ 2006; M endoncaj [200H |Nei> 
mayer et aLH2010"f Salamin et al.\ 2006 Son and Fisch, 



2005 ) for inertial confinement fusion ( Azechi et aZ.l |2006h 



based on th e high-energy density p lasma physics ( |Drake 



2009[ |2010[ |Norreys et~al\ |2009[ ). Plasma- like collec 



nanoscale magnetic switches. 

Collective interactions between an ensemble of degen- 
erate electrons and positrons/holes give rise to novel 
waves and structures in quantum plasmas. Studies of 
linear waves in a non-relativistic unmagnetized quantum 
plasma with degenerate electrons begun with the pio- 



tive behavior is well studied experimentally and theo- 
retically in solid state physics ( |Kittel[ [1996), in which 
metals and semi-conductors support both transverse op- 
tical modes, and longitudinal electrostatic modes, such as 
plasmons and phonons on electron and ion time-scales, 
and, in addition various lattice modes. Plasmons and 
phonons are usually probed by measuring the energy of 
electrons which have been passed through thin foils, or 
by laser scattering techniques. For example, the disper- 
sion relation of collective electron plasma waves has been 
measured for several metal specimen by using an elec- 



neering theoretical works of Bohm (19531; Bohm and 
Pines| (|1953|); |Klimontovich and Silin| ( |1952a|b[ |1961[ ) 



and Pines (1961), who studied the dispersion properties 



of high-frequency electron plasma oscillations (EPOs). 
The frequencies of the latter with an arbitrary elec- 



tron degeneracy have been found by Maafa (1993). In 



the theoretical description of the EPOs, Klimontovich- 
Silin and Bohm-Pines used the Wigner distribution func- 
tion ( Wigner 1932 ) and the density matrix approach to 



tron velocity analyzer of Mollenstedt type (Watanabe 



1956). Collective dispersive behavior of plasmons, in- 



cluding shifts in the plasmon frequency due to quantum 
effects, in solid-density plasmas have been observed by 



Glenzer et al. (2007) and Neumayer et al. (20101 using 



spectrally re s olved x-ray scatte ring techniques (Kritcher 
et al. 2008 Lee et al. 20091. In these experiments, 



demonstrate that in a quantum plasma with a Fcrmi- 
Dirac equilibrium distribution function for degenerate 
electrons, the frequency of the EPOs is significantly 
different from the Bohm-Gross frequency in a classical 
electron-ion plasma with non-degenerate electrons obey- 
ing the Maxwell-Boltzmann distribution function. The 
dispersion to the EPOs appears through the electron 
Fermi pressure and electron tunneling effects (|Gardner| 



and Ringhofer 1996 



Manfredi and Haas 



|Jungel et ^|2lX)6]|Manfredi[ |2005 



20011 Misra 



2009] |Shukla[ |2006 



3 



Shukla and Eliasson[ [20061 [20To| |Wilhelm[ [l97T| ). The 



quantum Bohm potential, responsible for electron tun- 
neling, appeared fir st in the quantum fluid description o f 
a sine le electron by IMadelungl (|1926|) and IBohm (11953 ) . 



For systems of degenerate electrons, different forms of 
the potential have been derived by using moments of the 



vestigating new aspects of collective interactions in dense 
quantum plasmas by means of non-relativistic quantum 
hydrodyna mic (|Gardner and Ringhofer 1996 Jungel et 
all [20061 |Manfredij |2005[ |Manfredi and Haas| |2001 



al 



Wigner equation (Ancona and Iafrate 1989 Iafrate et 
1981 ) and by using a variational approach QFeyn 



man and Kleinert 1986 Kleinert 1986 ) . They have been 



1998 



Shukl a and Eliasson[|2010[) and quantum kinetic (Bonitz 
|Kremp et aZ.111999] |Tsintsadze and Tsintsadze 



used i n qua ntum hydrodynamic (QHD) equations (|Wil- 



2009) equations. Models for non-ideal effects 
strongly coupled dense plasma ha v e been presented by 
Carruthers and Zachariasen (19831, Kremp et al. (20051, 



helm |1971 ) for modeling nano-devices (Ferry and Zhou 
1993 |Gardner and Ringhofer |1996[ ). More recently, La 



and Redmer and Ropke 
(WP) model ( |Hillery efal. 



(2010). The Wigner-Poisson 



rive a set of QHD equations (Manfredi 2005 Manfredi 



grangian approaches have been used to device efficient and Haas 2001 ) for ES waves in a quantum plasma. The 



computational algorithms for quantum systems (Loprc 



1984) has been used to de- 



ore and Wyatt 1999 Mayor et al. 1999). These and investigated by Haas et al. 



relation betwee n the QHD and kine tic models have been 

~ pOlObl. 



other methods for computational QHD using quantum 
trajectories have been nicely summarized in the textbook 
by |Wyatt| ( [2005l ). 

The quantum effects arc important for the dielectric 
and dispersive properties of a quantum plasma. The 
longitudinal and transverse dielectric constants of an 



isotropic quantum plasma were worked out by Lind 



hard ( 1954), Silin and Rukhadze ( 1961 ), and Kuzelev and 



Rukhadze ( 1999 1. Contributions of the electron spin and 
exchange interactions to the electromagnetic (EM) wave 
dispersion relations in an unmagnetized quantum plasma 



have been presented by Burt and Wahlquist (|1962) by 

:han- 



using a quantum kinetic theory. The quantum meet 
ical phase space distribution of Wigner ( 1932 ) has been 
further generalized by Brittin and Chappell ( 1962 ) for a 



system of charged particles including the quantized EM 
field and Green's functions involving correlations of dis- 
tribution functions and vector potentials. Kinetic models 
for spin-polariz e d plasmas have been devel oped by |Cow- 



ley et al. (|1986[),|Zhang and Balescu" (1988), and Balescu 
and Zhang| p988p . More recently, electron spin- 1/2 ef- 
fects in a quantum magnetoplasma have also been con- 



sidered by Brodin et al. (2008a) and Zamanian et al. 



(2010). The gauge problem in quantum kinetics has been 



treated by Stratonovich ( 1956 ) and Serimaa et al. ( 1986 1, 



which is important whenever the fields are not electro- 
static. In a quantum magnetoplasma, one finds that the 
external magnetic field significantly affects the dynamics 
of degenerate electrons, and that the thermodynamics 
and kinetics ( Steinberg 2000 1 in a quantum magneto- 



plasma are significantly different from those in an un- 



magnetized quantum plasma. Oberman and Ron ( 1963 ) 



derived the expression for the dielectric function for lon- 
gitudinal waves in a non-relativistic magnetized quantum 
plasma and discussed applications of their work to heav- 
ily doped semiconductors. Kelly ( 1964 ) studied the dis- 



persive properties of a magnetized quantum plasma by 
using the Wigner distribution function and the Maxwell 
equations. Finally, we mention that useful foundations 
for the theory of quantum plasmas are presented by [de] 
Groot and Suttorp ( 1972 ), while quantum kinetic models 



including the effects of spin are reviewed by Lee ( 1995 ) 



During the last decade, there has been a surge in in- 



ture (Manfredi and Haas 



2001 



The quantum na- 
Shukla and Eliasson[ 



2010) is manifested in the non-relativistic electron mo- 
mentum equation through the quantum statistical pres- 
sure, which requires the knowledge of the Wigner elec- 
tron distribution function for a quantum mixture of elec- 
tron wavefunctions, each characterized by an occupation 
probability. The quantum part of the electron pressure 



is also represented as a nonlinear quantum force ( Gard 



ner and Ringhofer 1996 Manfredi and Haas 2001 Wil 



helm[ |1971[ ) — V0s, where 4>b = —{h /2m eV / n^)V" % /n e 
is the Bohm potential, and m e and n e are the electron 
mass and electron number density, respectively.. Defining 
the effective wavefunction ij: — \J n e (r, t) exp[iS e (r, t) /K], 
where V5 e (r,t) = m e u e (r,t) and u e (r,t) is the elec- 
tron fluid velocity, the non-relativistic electron momen- 
tum equation can be cast into an effective nonlinear 
Schrodinger (NLS) equation (|Manfredi |2005 Manfredi 



and Haas[ [2001] [Shu kla , 2006] |Shukla and Eliasson[|2006[ 



2010), in which there appears a coupling between the 



electron wavefunction and the ES potential associated 
with the EPOs. The ES potential, in turn, is determined 
from Poisson's equation. One thus has the coupled NLS 
and Poisson equations, governing the dynamics of non- 
linearly interacting EPOs is a quantum plasma. Both 
non-relativistic QHD and NLS-Poisson equations exclude 
strong interactions among the quantum particles and 
electron exchange interactions ( Hohenberg and Kohn[ 
1964 Kohn and Sham 1965) between an electron and 
the background plasma particles (e.g. degenerate elec- 
trons and non-degenerate ions). However, it has turned 
out that the QHD and NLS-Poisson equations have been 
quite useful for studying linear and nonlinear plasma 



waves, as well as stability of quantum plasmas (Haas 
2005] [20071 l H aas et al.\ [2U031 |Manfredi[ [2005] |Manfredi| 



and HaasH200H|Shukla[[2006]|Shukla and Eliassonl|2006 



2010 ) at nanoscales involving the quantum force ( Gard- 
ner and Ringhofer[ 11996 Wilhelm 1971) and the quan- 



tum statistical pressure law for an unmagnetized quan- 
tum plasma with degenerate electrons. New effects also 
appear when one accounts for the potential energy of the 
electron spin-1/2 in a magnetic field (Brodin and Mark- 



hmd]|2007a|b|c[|Brodin et al 



Brodinj [2007] |Misra[ [2007] 



2008b, 2010; Marklund and 



2009 Misra and Samanta 



4 



2010 Shukla 



2007 2009 Takabayasi 1955). In fact, 



the QHD model for degenerate electrons in both non- 
relativistic dManfredil 120051 IManfredi and Haasl 12001 



Shukla[[2006HShukla and EliassonH2006| |2010[ ) anc 



rela- 

tivistic 1 Masood et al. |2010 1 quantum plasma regimes 
seems to provide an adequate description for probing 
some quantum collective interactions in compressed plas- 



mas 


(|Froula et al. 




2011 


Glenzer and Redmer| 2009 


Glenzer et al.\ |2007[ |Lee et al. 2009| |Neumayer et al. 


2010 


) due to the availability of ultrafast x-ray Thompson 



scattering spectroscopic techniques. 

In this Colloquium, we present the recent development 
of numerous nonlinear collective processes in a quantum 
plasma with degenerate electrons. We first describe the 
salient properties of quantum plasmas in which degen- 
erate electrons follow the Fermi-Dirac distribution. We 
then present the relevant equations for describing linear 
and nonlinear wave phenomena in quantum plasmas. Af- 
ter reviewing the linear properties of ES and EM waves, 
we proceed by presenting numerical results of the gov- 
erning nonlinear equations which reveal localization of 
ES and EM waves at nanoscales. Specifically, we discuss 
the formation and dynamics of nanostructures (e.g. ID 
quantum electron density cavity and 2D quantum vor- 
tices), as well as discuss the properties of 3D quantum 
electron fluid turbulence at nanoscales. Also presented 
are nonlinear interactions between intense EM waves and 
ESOs, which reveal stimulated scattering of EM waves off 
quantum plasma oscillations and trapping of light into a 
quantum electron density cavity. The effects of an ex- 
ternal magnetic field on linear and nonlinear wave phe- 
nomena in a quantum magnetoplasma are examined. Fi- 
nally, we highlight possible applications, as well as future 
perspectives and outlook of nonlinear quantum plasma 
physics. 



II. BASIC PROPERTIES OF QUANTUM PLASMAS 

Let us first summarize some of the basic properties of 
quantum plasmas that are quite distinct from classical 
plasmas. While classical plasmas are composed of non- 
degenerate plasma particles with low number densities 
and relatively high electron and ion temperatures, quan- 
tum plasmas have degenerate electrons and/or positrons 
with extremely high number densities and relatively low 
temperatures. The ions can usually be treated as non- 
degenerate plasma particles. Figure 1 depicts the plasma 
parameter regimes (the electron temperature versus the 
electron number density) under which quantum plasmas 
occur in different physical environments. 

Quantum mechanical effects start playing a significant 
role when the Wigner-Seitz radius (average inter-particle 
distance) a — (3 /47m) 1 / 3 is comparable to or smaller 
than the thermal de Broglie wavelength Xb = h/mVr, 
where m is the mass of the quantum particles (e.g. degen- 
erate electrons, degenerate positrons, degenerate holes), 
Vt = {ksT /m) 1 / 2 is the thermal speed of the quantum 
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particles, T is the temperature, m is the mass, and ks is 
the Boltzmann constant, i.e. when 



n\% > 1, 



(1) 



or, equivalently, when the temperature T is comparable 
to or lower than the Fermi temperature Tp = Ep/ks, 
where the Fermi energy is 



E F = — (37r 2 ) 2 /V/ 3 . 
2m 



(2) 



The relevant degeneracy parameter for the quantum 
plasma is 



T; 



f = \(^ 2 ) V3 (n\ 3 B ) 2/3 > 1. 



(3) 



For typical metallic densities of free electrons, n~5x 
10 22 cm~ 3 , we have Tp ~ 6 X 10 4 K, which should be 
compared with the usual temperature T. 

When the plasma particle temperature approaches Tp, 
one can show, by using a density matrix formalism 
( Bransden and Joachain 2000 1 , that the equilibrium dis- 



tribution function changes from the Maxwell-Boltzmann 
oc exp^—E/ksT) to the Fermi-Dirac (FD) distribution 
function 



Tfd — 2 



f— V 



1 + cxp 



E-fi 

knT 



(4) 



5 



where in the non-relativistic limit the energy is E = 
{m/2)v 2 = (m/2)(v 2 + v 2 + v 2 z ). The chemical poten- 
tial is denoted by fi. The parameter p/k B T is large and 
negative in the non-degenerate limit, and is large and 
positive in the completely degenerate limit. The equi- 
librium electron number density associated with the FD 
distribution function is 



n = / Tfd d 3 v — 



1 (2mk B T 
irh 2 



3/2 



Li 3/2 [-expfe)], 



( 5 ) 

where Li 3 / 2 is the poly- logarithm function, and £ M = 
fi/k B T. The completely degenerate limit corresponds 
to /i — » k B T F and T F 3> T. The relation between 



T F /T and is (Melrose 2008) -Li 3/2 [- exp(£ 



(4/3^)(T F /T) 3 / 2 . 

It us useful to define the quantum coupling parame- 
ters for electron-electron and ion-ion interactions. The 
electron-electron Coulomb coupling parameter is defined 
as the ratio between the electrostatic interaction energy 
Ei n t — e 2 /a e between electrons and the electron Fermi 
energy Ep e — k B T Fe , where e is the magnitude of the 
electron charge and a e = (3/47rn e ) 1 / 3 is the mean inter- 
electron distance. We have 



E; 



E Fe 



0.3 



1 



2/3 



0.3 



k B T Fc 



where X Fe = V Fe /uj pe , V Fe = (2E Fe /m e ) 1/2 = 
(h/m e )(3ir 2 n e ) 1 / 3 is the electron Fermi speed, and uj pe — 
(47rne 2 /m e ) 1 / 2 the electron plasma frequency. Further- 
more, the ion-ion Coulomb coupling parameter is Ti = 
Zfe 2 / diksTi, where is the ion charge state, dj = 
(3/47rni) 1 / 3 the mean inter-ion distance, and Tj the ion 
temperature. 

Since Y e for metallic plasmas could be larger than 
unity, it is of interest to enquire the role of inter-particle 
collisions on collective processes in a quantum plasma. It 
turns that the Pauli blocking reduces the collision rate for 
most practical purposes ( |ManfredI 2005 Son and Fisch 



2005 ) . Due to Pauli blocking, only electrons with a shell 



of thickness k B T about the Fermi surface suffer colli- 
sions. For these electrons, the electron-electron collision 
frequency is proportional to k B T '/%. The average colli- 
sion frequency among all electrons turns out to be ( Man- 
rredil[2005| 



k B T 2 



Fe 



(7) 



Typically v ee <C w pe when T < T Fe , which is relevant 
for metallic electrons. On the other hand, the typical 
timescale for electron-ion (lattice) collisions is r ei ~ 10 
fs, which is one order of magnitude greater than the elec- 
tron plasma period. Accordingly, a collisionless quantum 
plasma regime is relevant for phenomena appearing on 
the timescale of the order of a femtosecond in a metallic 
plasma. 



In compact astrophysical objects such as white dwarf 
stars, the mean distance n e between electrons become 
comparable to the Compton length Xq = h/m e c, and 
accordingly the speed of an electron on the Fermi sur- 
face becomes comparable to the speed of light c in vac- 
uum, so that one has to take relativistic effects into ac- 
count. Relativistic degenerate electrons are found in the 
core of massive white dwarf stars (Koester and G. Chan- 
mugam 1990 Shapiro 1983[ ), aptly named due to their 
very low luminosities yet high surface emissivities, which 
are compact bodies with radii < 10 -2 Rq and masses 
typically < M . Consequently, the average electron 
number densities are quite high (~ 10 30 cm -3 ). Since 
electrons are Fermions, only one electron can occupy a 
given quantum state (position, spin) . In a simplified pic- 
ture, each electron will on average occupy a vol ume l/n e . 
Then, by Hei s enber g's uncertainty principle (Bransden 
and Joachain 2000), AxAp < 7i/2, the mean momen- 

1/3 

turn of electrons can be estimated to be p x « fin e ■ If 
electrons are non-relativistic, the velocity of the electron 
is ~ Px/m-e — hnl^ 3 /m e ; however, if electrons are rela- 
tivistic, their velocity will be close to c. Now the electron 
pressure, as it is for a simple gas, is the momentum trans- 
fer per unit area, or P e = (momentum) x (velocity) x 
(number densit y) . For non-relativistic electrons, we have 



(6) (iGurskyl [1976J) P e = hnl /3 (hnl /3 /m e )n e = K z n° e /ci /m, 



^2,5/3, 



On the other hand, when electrons are relativistic, the 

4/3 



relativistic electron pressure is P er = Hn e cn e = Kcn e 
In the past, |Chandrasekhar ( 1931a|b 1935, 1939) and 
others presented a rigorous derivation of the electron 
pressure Pc for arbitrary relativistic electron degeneracy 
pressure in dense matter. It reads 



Pc 



3h 3 



mjc 5 



/(&), 



(8) 



where /fe) = £ c (2£ c 2 - 3)(1 + £?) 1/2 + Ssinh -1 ^,), & = 
p c e/m e c, and p c = (3/i 3 n e /87r) 1 / 3 is the momentum of 
an electron on the Fermi surface. In the non-relativistic 



limit £ c <C 1, we have (Chandrasekhar 1935 1939) 



P„ 



2/3 



5m e 



h 2 nl' 3 . 



(9) 



while in the ultra-relativistic limit £ c ^> 1, the degenerate 
electron pressure reads (Chandrasekhar 1931a) 



P„ 



.1/3 



-heni' 3 = -hent' 3 . 



(10) 



Thus, the intuitively obtained formulas of Gursky ( 1976 ) 



for non-relativistic and ultra-relativistic pressures for de- 
generate electrons are in agreement with those deduced 
from the pressure formula (8) for an arbitrary relativistic 
electron degeneracy pressure. 

In his Nobel Prize winning papers on the structure 
of compact stars, Chandrasekhar ( 1931a|b I balanced 



the gradient of the ultra-relativistic electron degener- 
acy pressure P u /R and the gravitational force G — 



6 



(GM/ R 2 )n e m n , where G is the gravitational constant, 
M and R are the mass and radius of a star, respectively, 
m n the mass of the nuclei (n e m n = M/R 3 ), to deduce the 

critical mass of a star M c — (hc/G) 3 ^ 2 m~ 2 w 1.4 M Q , 
where M is the solar mass. The interior of white dwarf 
stars usually consists of fully ionized helium, carbon, and 
oxygen, which approximately consist of equal amounts of 
protons and neutrons. Hence, the effective mass of the 
nuclei can be taken to be the proton mass plus the neu- 
tron mass. Since M c is independent of density, it means 
that this mass is obtained independent of radius. This 
is the limiting mass; more massive stars cannot be sup- 
ported by electron degeneracy pressure no matter how 
small they are. This was the discovery of Chandrasekhar; 
that the pressure dependence on density changed in go- 
ing from nonrelativistic to relativistic conditions and, as 
a consequence, there arose a finite limit to the mass of a 
star with ultra-rclativistic degenerate electrons. 



III. MODEL EQUATIONS FOR QUANTUM PLASMAS 

In quantum systems, the Dirac and Maxwell equations 
are often used to study the dynamics of a relativistic 
quantum particle/Fermion (electrons and positrons) in 
the presence of intense electromagnetic fields. Quan- 
tum particles have spin. For example, an electron spin 
s = 1/2 is an intrinsic property of electrons which have 
an intrinsic angular momentum characterized by quan- 
tum number 1/2, and a magnetic moment for individual 
electrons. In fact, the relativistic Dirac equation provides 
a description of quantum particles (with spin) under the 
action of the electromagnetic fields. The spin of elec- 
trons (and positrons)-which have the spi n- 1/2 has be en 
introduced through Dirac's Hamiltonian (Diracl 1981) 



where 



H a — 



h —n iq a h , . _ _ . . s 

V 2 - (A- V + V- A) 



2m a 
2m a c 2 



2m a c 
q a (j) - fj, a ■ B, 



(13) 



is the Hamiltonian, and i/j a (r,t,a) is the wavefunction 
of the single quantum particle species a with the spin 
s = (l/2)er (|<t| = 1), and q a = —e (+e) for electrons 
(positrons). The last term in (13) is the potential en- 
ergy of the magnetic dipole in the external magnetic field, 
the magnetic moment of which is /i Q = (q a h/2m a c)(T = 
/j,b (T , where fis — q a h/2m e c is the Bohr-Pauli magne- 
ton and <r the spin-operator of a single quantum particle 
( Landau and Lifshitz[ |1998a[ ) 



By using the Madelung representation (Madelung 



1926) for the complex wavefunction ip a , viz. 



ip a (r, t, a) = * Q (r, t, a) exp (iS a /K) . 



(14) 



where ^ a (r,t,a) and S a (r,t,a) are real, in the Pauli 
equation ( 12), we obtain the quantum Ma delung fluid 
equations ( Tsintsadze and Tsintsadze 2009 ) 



V • (n a p a /m a ) = 0, 



and 



dt 



Uq, x B\ 
q a |E+^ +Fq+F s , 



where we have denoted 



2m a 



V 2 



(15) 



(16) 



(17) 



and 



Tl = ca s ■ (p e + -A) — ecf) + f3m e c 2 , 



(11) 



where p e = —iTiV is the momentum operator, and ct s 
and (3 are the Dirac matrices. The three Cartesian com- 
ponents oij (j — 1,2,3) of a s are usually cons tructed 
with help of the Pauli spin matrices a x , a y and a z ( Brans- 
den and Joachain 2000 1 . The corresponding wave func- 
tions ip are four-component spinors. The magnetic field 
is B = V x A, where A and <j> are the vector and scalar 
potentials, which are determined from the Maxwell equa- 
tions. 

In the non-relativistic limit, the Pauli equation 



( Berestetskii et al. 1999) in the presence of the electro- 



magnetic fields describes the dynamics of a single quan- 
tum particle. It reads (Tsintsadze and Tsintsadze 2009) 



ift 



dt 



H a i/j a 



(12) 



F s = AisV (a • B) 



(18) 



Here n a = l^a) 2 is the probability density of finding a 
single quantum particle with a spin s at some point in 
space, p Q = V5 Q — q a A/c is the momentum operator 
of a quantum particle, d/dt = (d/dt) + u a ■ V, u Q is the 
velocity of a quantum particle, and E = — V0— c~ 1 dA/dt 
and B = V x A. 

The spin force F s in a quantum magnctoplasma can 
also be written as ( Brodin and Marklund"| 2007a|b|c 
Brodin et al] |2008bl Marklund and Brodinl [20071) 



F s=AlB tanhf^-) VB, 
\k B T a J 



(19) 



where B = |B| and tanh(£) = B 1/2 (^), with the Bril- 
louin function with argument "1/2" describing particles 
of spin-1/2. The Langevin parameter tanh(£) accounts 
for the macroscopic magnetization of electrons due to 
the electron thermal agitation and electron-electron col- 
lisions. 
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A. The Schrodinger and Wigner-Poisson Equations 

The quantum A-body problem is governed by the 
Schrodinger equation for the A-particle wavefunction 
tp(qi, 92, • • • , 3jv)j where qj — (r^, Sj) is the coordinate 
(space, spin) of the particle j, each particle associated 
with energy £j . A drastic simplification occurs if one ne- 
glects the correlation between the particles at every order 
in Tq and describes the full wavefunction as the product 
of the single particle wavefunctions. For identical quan- 
tum particles, the A-p article wavefunction is given b y 
the Slater determinant ( Bransden and Joachain 2000 ) 



ip(qi, 92, ■ ■ • , 9jv) 

i>i(qi,t) ip2(qi) ■■■ -iM<7i) 
4>i(q2,t) -02(92) ••• V'ivte) 



7 A! 



4>l{qN) ^2(9^) ••• ^7v(9Af) 



(20) 



which is antisymmetric with respect to an interchange 
of any two particle coordinates. This property is re- 
quired by the Pauli exclusion principle under the sec- 
ond quantization procedure for a system of A identical 
non-relativist ic quantum particles. Accordingly, tp van- 
ishes if two rows are identical, i.e. two identical quan- 
tum particles cannot occupy the same state. Example 

(A = 2): V(<7i,<72) = ^[M<h)Mto) - M&)Mqi)] 

so that tp(q2,qi) = —"0(91,92) and ip(qi,qi) — 0. In the 
zero temperature limit, all energy states up to the Fermi 
energy level are occupied, while no energy states above 
the Fermi level are occupied. 



To c aptur e coll ective effects in quantum plasmas, Haas 
et al. (2000) and Anderson et al. (20021 used the time- 



dependent Hartree model where electrons are described 
by a statistical mixture of A pure states, where each 
wavefunct ion ?/>,-, j = 1, . . . , A obeys the Schrodinger 



equation (Anderson et al. 2002 1 



^ + 2^ V ^ +e ^- ' (21) 
which is coupled with Poisson's equation 



Aire 




ZiTli 



(22) 



of like-spin to avoid each other. Each electron of a 
given spin is consequently surrounded by an "exchange 
hole", a small volume around the electron which like- 
spin electrons avoid. For the study of magnetic ordering 
in quantum dots (QDs) doped with magnetic impurities, 
Eqs. (21) and (22) must also be enlarged by including 
a 3D QD confining potential and a Vosko-Wilk-Nusair 



spin dependent exchange-correlation potential ( Dharma- 



wardana and Perrot 1995). Hence, the self-consistent 



model will go far beyond the Kohn-Sham's description 



(Kohn and Sham 1965) for treating the dynamics of 



correlated electrons in electron clusters, accounting for 
electron-exchange and electron-correlation effects. In the 
presence of time-dependent potentials, the properties and 
dynamics of many-electron systems can be investigated 
by using a ti me-dependent functional theory ( |Runge and| 



Gross 1984) 



However, in a non-relativistic quantum plasma with 
an ensemble of degenerate electrons, it is more appropri- 
ate to use the quantum stati stical theory in volving the 
Wigner distribution function (Wigner, 1932) 



x J exp(im e v • R//i)V>*(r + R/2)0(r - R/2) d 3 R, 

(23) 

where the asterisk denotes the complex conjugate. Eq ua- 
tion (23) has also been used by Moyal (Moyal 19491 for 
studying the dynamics of electrons in a quantum system. 

For electrostatic interactions in a quantum plasma, the 
Wigner-Poisson equations, to a leading order (in the limit 
of weak quantum coupling parameter T e ), can be written 
as 



df w 
dt 



iem% 

V • V^ = -(2^' , ' 

R \ ,/ R 

x + — ,t — <p x — — ,t 



m s (v-v')-R/h 

/ M (x,v',t) d 3 Rd 3 v' 
(24) 



and 



X7 z (h = 47re 



f w d 3 v - Z t n, 



(25) 



where rii is the ion number density (to be obtained from 
the hydrodynamic equations for non-degenerate ions, to 
be discussed later) , and the electrostatic potential aris- 
ing from the charge distribution of A electrons. Equa- 
tions (21) and ( 22) have been u s ed to study stream- 
ing instabilities (Anderson et al. 2002) and other ki- 



netic effects in a quantum system composed of an en- 
semble of electrons. Within the Hartree-Fock model, 
Eq. (21) can be further generalized by including the 
electron-exchange term resulting from the Pauli exclu- 
sion principle. The effect of exchange is for electrons 



B. The QHD Equations 



The non-relativistic QHD equations (Wilhelm 1971) 
have been developed in co ndensed matter physics (Gard- 
ner and Ringh ofer] 1996 ) and in plasma physics ( Man- 
fredil 120051 IManfredi and Haasl 120011). The non- 



relativistic QHD equations are composed of the electron 
continuity equation 



dn e 
~dt 



V ■ (n e u e ) = 0, 



(26) 
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the electron momentum equation ( Wilhelm 1971 1 

m e ^ + u e .Vu e J =eV0- ^-VP e +F Q , (27) 

and Poisson's equation 

V 2 = 47re(n e — ZiUi). 



(28) 



In a quantum plasma with non-relativistic degenerate 
electrons, the quantum statistical pressure in the zero 



electron temperature limit can be modeled as ( Crouseilles 
et al.\ 120081 IManfredi and Haasl 120011) 



P, 



T/2 / \ (D+2)/D 



n 



(29) 



where D is the number of space dimension of the system, 
and V Fe — {fr/m e )(3iT 2 n e ) 1 / 3 the electron Fermi speed. 



C. The NLS-Poisson Equations 

For investigating nonlinear properties of dense quan- 
tum plasmas, it is appropriate to work with a NLS equa- 
tion. Hence, by introducing the wavefunction 



ip(r,t) = Vn e (r,i) 



exp 



.S e (r,t) 



(30) 



where S e is defined according to m e u e =X7S e and n e — 
\ip\ 2 it c an be shown that (27) can be cast into a N LS 
equation (Manfredi 2005 Manfredi and Haas 2001) 



ih 



dip h 2 

— H 

dt 2m, 



-v v + e-H> 



m eVpe 

2n 2 



|^| 4/ ^ = 0, (31) 



where the electrostatic field <p is determined from Pois- 
son's equation 



V 2 = 47re(|V 2 | - ZiO. 



(32) 



We note that the third and fourth terms in the left-hand 



side of Eq. (31) represent the nonlinearities associated 
with the nonlinear coupling between the electrostatic po- 
tential and the electron wavefunction and the nonlinear 
quantum statistical pressure, respectively. 



where k is the wavenumber and uj pe = (47rrioe 2 /m e ) 1 / 2 is 
the electron plasma frequency. Here the ions are assumed 
to be stationary. 

One can identify two distinct dispersion effects from 
(33): One long wavelength regime with V Fe ^ Tik/2m e , 
and the other short wavelength regime with V Fe < 
hk/2m e . These two regimes are separated by the crit- 
ical wavenumber 



2tt 



v crit 



m e V F e 



-1/3 



(34) 



It should be mentioned here that the quantum dispersion 
effects associated with the EPOs have recently been ob- 



served in a compressed plasma (|Froula et al. 


pOll Glen- 


zer et al. 


2007 


Neumayer et al. 


2010). In compressed 



plasma experiments, powerful x-ray sources are employed 
for accessing narrow bandwidth electron plasma wave 
spectral lines via collective Thomson scattering in which 
powerful light scatters off electron density fluctuations. 
We note that the dispersion relation for EPOs in the fi- 



nite electron temperature limit is given by ( Thiele et al. 
20081 



cu 2 = lu 2 + 3/c 2 y| e (l + 0.088n e A 3 ) + 



h 2 k 4 

Am 2 



(35) 



where Vr e — (kBTe/me) 1 ^ 2 is the electron thermal speed 
and A e = ^/2~TTh/y/m e kBT e the thermal (De Broglie) 
wavelength. 

As mentioned in the Introduction, in the past many 
authors derived the dielectric constant for the high- 
frequency (in comparison with the ion plasma fre- 



quency) ES waves i 


Bohm[ 


1953 


Bohm and Pines 


1953 


Klimontovich and £ 


ilin 1952a|b| 


Lifshitz and Pitaevskii 



1981[) and the refractive index for EM waves (Burt and 



Wahlquist 1962) by using a quantum kinetic theory 
based on the Wigner and Poisson-Maxwell equations in a 
quantum plasma. In the following, we briefly discuss the 



well known results for the ES (|Bohm 


1953 




Bohm and 


Pines 1953 Klimontovich and Silin 




1952a 


b 


) and EM 


(Burt and Wahlquist 1962) waves in an unmagnetized 



quantum plasma. 

The dielectric constant for ES waves in a plasma 
with completely degenerate electrons reads (Lifshitz and 
Pitaevskii 



1981) 



IV. LINEAR WAVES IN QUANTUM PLASMAS 

A. Electron Plasma Oscillations (EPOs) 

Linearization of the NLS-Poisson Equations (31) and 
(32) around the equilibrium state and combining the re- 
sultant equations, we obtain the frequency ui of the EPOs 
( Bohm[ [l953[ [Bohm and Pines 1953; Kli montovich and] 
Silin[ |1952a|bt 



h 2 k^ 1/2 



Am 2 



(33) 



3w 2 



2k 2 V 2 c 



[l-g(u+)+g(u-)], (36) 



where w± = u> ± fik 2 /2m e , and 



9(u±) = 



k 2 V 2 Fe 



2hkV Fe 



-log 



u)± + kVi 



Ft: 



oj± — kVp 



(37) 



Assuming that the phase velocity (ui / k) of the ES wave is 
much larger than V Fe , we obtain by setting D e (w,k) = 
the frequency of the EPOs, given by (36). On the 
other hand, in the semi-classical limit, viz. ?i|k| <C p Fe = 
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h(3ir 2 n e ) 1 ^ 3 , we have (Lifshitz and Pitaevskii 1981 1 from 
Eq. (36) 



D e (w,k) = l 



1 - 



UJ 

2kV Ft 



log 



uj + kVp t 

UJ — kVp t 



which in the short wavelength limit, viz. kVp e 3> 
yields the so-called an electron thermal quasi-mode 



pe j 



( |Klimontovich and Silin] |1952a|b| [l961[ ) 

uj = kV Fe [1 + 2exp (-2fc 2 A 2 - 2)] , 

^Fe/VS is the Thomas- Fermi screening 



(39) 



where A s ~ 
length. 

Furthermore, when uj = 0, the expression (38) as a 
function of k has a Kohn singularity at hk = 2pj? e = the 
diameter of the Fermi sphere. Here we have 



D e (0,k) = l 



2whE f 



[i-eiog(i/|CD], (40) 



where £ = (hk — 2pp e )/2pp e and |£| < 1. In a quan- 
tum plasma, with D(0,k) given by (38), the potential 
distribution (p(r) around a stationary test charge q t is 



4:irq t f exp(ik • r)d 3 k 



^ (2ir) 3 J & 2 A(0,k) 
which gives ( |Else et al.\ |2010[ ) 

cos(2fc^r) 



ip(r) 



Qt 



(2 + 3/7 2 



(41) 



(42) 



- PFe/h. We note that 
3 cos(2fci?r), consider- 



where r\ — huj pe / 'AksTpe and kp 
(42), which is proportional to r 
ably differs from the Debye-Hiickel shielding potential 
that is proportional to r" 1 exp(— r/\p, e ) in a classical 
plasma with the Maxwell-Boltzmann electron distribu- 
tion function. Here Xp> e is the electron Debye radius. We 
further note that the shielding of a moving test charge in 
an unmagnetized quantum plasma has been investigated 
by Else et al. (20101 both analytically and numerically. 



B. Ion Plasma Oscillations (IPOs) 



where is the ion number density, the ion fluid veloc- 
ity, mi the ion mass, pi = riiVtii the ion mass density, 
the adiabatic index for the ion fluid, r m the viscoelastic 
relaxation time for ions, rj and £ the bulk ion viscosities. 
The viscoelastic equation (44) for strongly systems has 
been successfully used (Ichmaru and Tanaka 1986 Kaw 



and Sen 1998) for investigating collective processes in 
classical plasmas with non-degenerate plasma particles. 

The ions are coupled with degenerate electrons by the 
space charge electric field E = — V0. For low-phase ve- 
locity (in comparison with the electron Fermi speed) ES 
waves, we can neglect the inertia of the electrons to ob- 
tain 



9 ri 



Ti 2 n e 



n e V0- \7n 5 J 3 , 

5 m„ Am 



V 2 n, 



0, (45) 



for a quantum plasma with weakly relativistic degener- 
ate electrons, while for a quantum plasma with ultra- 
relativistic degenerate electrons, we have 



n e V0- -hcVnj /3 = 0. 



(46) 



Due to the ion inertia, one has new dielectric con- 
stants for the low-frequency (in comparison with the elec- 



tron plasma frequency) ES waves (Eliasson and Shukla 


2008a] pushtaq and Melrose[ |2009 Pines 


1963 Pines| 


and Nozieres 




1989| |Shukla and Eliasson 




2008b 1. In a 



quantum plasma with non-relativistic degenerate elec- 
trons with uj 2 <C k 2 Vp e + ?i 2 fc 4 /4m 2 , we can linearize 
(28), (43), (44), and (45), Fourier transform them, and 
combine the resultant equations to obtain 



A(w,k) = 1 



3uj 2 e 



UJ 



k 2 V 2 e 



h 2 k A /Am 2 e 



pi 



(47) 



where uj pi = (Aim^Zfe 2 /m,) 1 / 2 is the ion plasma fre- 
quency, and Of = uj 2 — r jik 2 V^ i + iujk 2 i]*/(l — iujT m )), 
with Vt% — (kBTi/nii) 1 / 2 and ?y* = (£ + 4rj / 3) / m^rto . 
On the other hand, in a quantum plasma with ultra- 
relativistic degenerate electrons, we have from (28, (43), 
(44), and (46) 



We now focus our attention on the effect of the dynam- 
ics of non-relativistic and non-degenerate ions in an un- 
magnetized quantum plasma. The dynamics of strongly 
coupled ions is governed by the ion hydrodynamic equa- 
tions composed of Poisson's equation (28), and the con- 
tinuity and momentum equations. The latter are 

' ^ • (n^) = 0, (43) 



and 



dt 
d 

8t + 



Zi.e. 



'-V<j>- 



V(V • ik) = 0, 



(44) 



A(w,k) = 1 



k 2 C 2 



(48) 



where we have denoted C 2 t 



c^Xcnj/ 3 , and Ac = h/m e c 
is the Compton length. By setting Djfu^k) = 0, we 
obtain the frequencies of the IPOs. For the case with 
non-relativistic degenerate electrons we have 



uj 2 = 7ik 2 V Tl + 



k 2 r\^ 



uj pi rv /\ Th 

(i + k 2 x 2 Th y 



(49) 



while for the case with ultra-relativistic degenerate elec- 
trons the result is 



uj 2 = ll k 2 V 2 l + ^± 



u pi 2k 2 X 2 h 

(i + k 2 x 2 y 



(50) 
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where we have assumed ujT m <C 1 and denoted Xth = 

1/2 

[A 2 + h 2 k 2 / 4mlujp e ] , and X h = Cn/uj pe . The domain 
of validity of the hydrodynamic description for the ions 
in the context of ion oscillations in a weakly relativistic 
den se pla sma has also been recently discussed |Mithen ei| 

oDpoul). 



Melrose and Mushtaq ( 2009 ) and Mushtaq and Melrose 



(2009) have presented Landau damping rates for both 



electron and ion plasma waves in an unmagnetized dense 
quantum plasma. The imaginary parts of the dielectric 
constants can be used to calculate the structural form 



generate electrons. 



factor ( Ichimaru . 1982 ) in quantum a plasma with de- 



Shukla and Eliasson (2008b) used the dielectric con- 



stant (47) without the quantum statistical pressure term 
(viz. the V^ e -term) to investigate the screening and 
wake potentials around a test charge in an electron-ion 
quantum plasma. They found a new screening potential 
(iShukla and Eliassonl |2008b|) 



— exp(— k„r) cos(fc„r), 
r 



(51) 



and the wake potential 



qt 



"0*| 



cos 



(z - U t) 

u 



where k q = \/2/ y/h/m e uj pe is the quantum wave num- 
ber, and r = [x 2 + y 2 + (z — UQt) 2 ] 1 / 2 the distance from 
the test charge moving with the speed Uq along the z 
axis in a Cartesian co-ordinate system. The wake po- 
tential (52) behind a test charge arises due to collective 
interactions between a test charge and the ion oscilla- 
tion with the frequency uj k w uj pi kj_/(k 2 _ + k 2 ) 1 / 2 , with 

k z < k qi kj_ = (fc 2 + k 2 ) 1 / 2 . We note that the Shukla- 
Eliasson (SE) exponential cosine-screened Coulomb po- 



tential 



has a minimum of 



-0.02q t k q at 



3k~ , similar to the Lennard- Jones potential for atoms. 
The SE screening potential ip se , which is independent of 
the test charge speed uo, is different from the Yukawa 
screening potential (?t/V) exp(— r/X s ) that is valid in the 
limit Vp e ^> hk/2m e . Recentl y , sever al authors (|Ghoshal 
" 2009a|b 



and Ho 



C. High- Frequency EM Waves 

Finally, we turn our attention to the high-frequency 
(HF) EM waves in an unmagnetized quantum plasma. 
Noting that HF-EM waves in the latter do not give rise 
to any density perturbations, we have the EM wave fre- 
quency 



1 1 2 2 i 2 \l/2 



(54) 



However, consideration of the electron spin current and 
electron exchange potential contributions in a quantum 
plasma gives rise to additi onal contributions to the r e- 
fractive index N. We have (Burt and Wahlquist 1962) 



= Nml- 



pe 



u 2 pe h 2 k 2 



1 



-Ki 



1 



, (55) 



which includes the electron spin correction, and is valid 
at zero temperature. Here TiKp = (2m e Ep e ) 1 ^ 2 is the 
momentum of degenerate electrons at the Fermi surface, 
the (l/5)Kp term is related to the leading quantum term 
from the ordinary transverse current, and the k 2 /4 term 
arises from the electron spin interactions. On the other 
hand, the EM wave dispersion relation, which accounts 
for the electron exchange potential and discards the spin 



(52) correction, reads (Burt and Wahlquist, 1962) 



2^,2 



k 2 c 



N w 1 



u 2 pe h 2 k 2 K 2 
5m 2 w 4 



3^ 2 k 2 r 
40^ ^ 



V. QUANTUM DARK SOLITONS AND VORTICES 

Let us now discuss nonlinear properties and dynamics 
of ID quantum dark solitons (characterized by the lo- 
cal electron density depletion associated with a positive 
potential) and 2D azimuthally symmetric electron vor- 
tices in an unmagnetized quantum plasma (Shukla and 
Eliassonl 2006) with immobile ions. The assumption of 



stationary ions is justified because we are looking for the 
nonlinear phenomena on a timescale much shorter than 
the ion plasma period. 



Xia et al.\ [2010| have used the SE |2006| |Shukla and Eliasson[ |2006[ ) 



We use the normalized NLS-Poisson equations ( Shukla 



potential to study doubly excited resonance states of 
Helium and hydrogen atoms emb edded in a quantum 
plasma ( Ghoshal and Hoj 2009a|b I, and lattice waves in 
2D hexagonal quantum plasma crystals (Xia et al. 2010 1. 



and 



Furthermore, by using D t from (48), one can deduce 
potential distributions around a moving test charge in 
a quantum plasma with ultra-relativistic electrons. We 
have 



^IV 2 * + (pV - \^\ 4/D y = 0, 



vVH*! 2 -!, 



(57) 



(58) 



(r, z) = — exp 
r 



(It 



A c J \z — uo*| 



(z - u t) 



where Ac = C^/w pe and L c 
M = u /C h . 



X r M 2 -1 



1/2 



(53) 
> 0, with 



where the time and space variables are in units of 
h/ksTpe and the electron Fermi-Thomas screening 
length Xtf, respectively. Furthermore, we have denoted 

VP = ip/y/no, if — etp/kBTpei and A = 2irn^ ? ' 'e 2 /ksTpe- 
The system ( 57 ) and ( 58 ) is supplemented by 



dt 



(59) 
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where = — Vip. Equ ations (f5T|) — ([59]) ha ve the follow- 
ing conserved integr als ( |Shai kh and Shukla 2007 Shukla 
and Eliasson 2006): the number of electrons 



N = J \V\d 3 x, 
the electron momentum 

p = ~iJ #*Wd 3 :r, 
the electron angular momentum 

L = -i J **r x W<i 3 x, 
and the total energy 



£ 



-_4**V 2 * 



D 



(2 + D) 



1^1(2+4/13) 



(60) 



(61) 



(62) 



d 3 x. 

(63) 




FIG. 2 The electron density \W\ 2 (the upper panel) and ES 
potential tp (the lower panel) associated with a dark soliton 
supported by the system of equations ( |64[ | and ( |65[ ) , for A = 5 
(solid lines), A = 1 (dashed lines), and .4 = 0.2 (dash-dotted 
line). After |Shukla and Eliasson| (|2006|). 



A. Quantum Electron Cavity 

For quasi-stationary, ID nonlinear structures moving 
with a constant speed vq, one can find solitary wave so- 
lutions of Eqs. (57) and (58) by introducing the ansatz 
"J = W(£) exp(iK s x — iCl s t), where W is a complex- 
valued function of the argument £ = x — vot, and K s 
and fl s are a constant wavenumber and frequency shift, 
respectively. By the choice K s — vq/2A, the coupled 
system of equations (57) and (58) can then be written as 



d 2 W 



and 



d£? 



tpW 
~A~ 



A 



0, 



1, 



(64) 



(65) 



where A = (£l s /A) — v^/AA 2 is an eigenvalue of the sys- 
tem. From the boundary conditions |W| = 1 and tp = 
at |£| = oo, we determine A = 1/A and fl s = 1 + Vq/4A. 
The system of Eqs. (64) and (65) admits a first integral 
in the form 



H h =A 



dW 


2 


1 (d v \ 






2 \di) 



2 

3 ~ 



(66) 



1 and tp = at 



where the boundary conditions \W\ 
|£| = oo have been employed. 

Figure [2] shows profiles of \W\ 2 and tp obtained numer- 
ically from ( [64] ) and (65) for a few values of A, where 
W was set to —1 on the eft boundary and to +1 on the 
right boundary, i.e. the phase shift is 180 degrees be- 
tween the two boundaries. The solutions are in the form 




FIG. 3 The time-development of the electron density |^| 2 
(left-hand panel) and ES potential tp (the right-hand panel), 
obtained from a simulation of the system of equations 
(57|) and (58]). The initial condition is * = 0.18 + 
tanh[20sin(x 7l0)]exp(i J g s x), with K a = v /2A, A = 5 and 
v = 5. After iShukla and Eliassonl (12006 1. 



of dark solitons, with a localized depletion of the electron 
density N e = \W\ 2 , associated with a localized positive 
potential. Larger values of the quantum coupling pa- 
rameter A give rise to larger-amplitude and wider dark 
solitons. The solitons localized "shoulders" on both sides 
of the density depletion. 

A numerical solution of the time-dependent system of 
Eqs. ( 57 ) and ( [58| is displayed in Fig. |3j with initial con- 
ditions close (but not equal) to the ones in Fig. [2] Two 
very clear and long-lived dark solitons are visible, asso- 
ciated with a positive potential of tp ~ 3, in agreement 
with the quasi-stationary solution of Fig. [2] for A = 5. 
In addition there are oscillations and wave turbulence in 
the time-dependent solution presented in Fig. [3] Hence, 
the dark solitons seem to be robust structures that can 
withstand perturbations and turbulence during a consid- 
erable time. 
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FIG. 4 The electron density l*^ 2 (upper panel) and ES po- 
tential ip (lower panel) associated with 2D electron vortices 
supported by the system ( 67 1 and ( 68 1 , for the charge states 
s = 1 (solid lines), s = 2 (dashed lines) and s = 3 (dash- 
dotted lines), with A = 5 in all cases. After Shukla and 
Eliasson|(|2006l. 



B. Quantum Electron Vortices 

For two-dimensional (D = 2) EPOs in quantum plas- 
mas, one can look for quantum vortex structures of the 
form \Jj = ip(r) exp(is9 — ifl v t) , where r and 6 are the po- 
lar coordinates defined via x = rcos(6) and y = rsin(6'), 
Q v is a constant frequency shift, and s — 0, ±1, ±2, . . . 
for different excited states (charge states). With this 
ansatz, Eqs. (571 and (58) can be written as, respec- 
tively, 



dr 2 



and 



1 d 

r dr 

cP_ 
dr 2 



1 d 

r dr 



^ = 0, (67) 



(68) 



where the boundary conditions ip — 1 and ip — dip /dr = 
at r = oo determine the constant frequency Q v = 1. Dif- 
ferent signs of the charge state s describe different rota- 
tion directions of the quantum vortex. For s / 0, one 
must have ip = at r = 0, and from symmetry consid- 
erations one has d(p/dr = at r = 0. Figure [4] depicts 
numerical solutions of Eqs. (67) and (68) for different 
values of s and for A = 5. Here a quantum vortex is char- 
acterized by a complete depletion of the electron density 
at the core of the vortex, and is associated with a positive 
ES potential. 

A time-dependent solution of Eqs. ( [57] ) and (58) in 
two-space dimensions for singly chargecl[s = iljclec- 
tron vortices is shown in Fig. [5j where, in the initial con- 
dition, four vortex-like structures were placed at some 
distance from each other. The initial conditions were 
such that the vortices are organized in two vortex pairs, 
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FIG. 5 The electron density |^| 2 (left panel) and an arrow 
plot of the electron current i (WVW* — \E r * V^) (right panel) 
associated with singly charged (|sj — 1) 2D electron vortices, 
obtained from a simulation of the time-dependent system of 
equations (57) and (58 1, at times t — 0, t — 3.3, t = 6.6 



and t — 9.9 (upper to lower panels), with A = 5. The singly 
charged vortices form pairs and keep their identities. After 
Shukla and Eliassonl (12006 1. 



with si = +1, s 2 = —1, S3 = —1, and s 4 = +1, seen 
in the upper panels of Fig. [5j The vortices in the pairs 
have opposite polarity on the electron fluid rotation, as 
seen in the upper right panel of Fig. [5j Interestingly, 
the "partners" in the vortex pairs attract each other and 
propagate together with a constant velocity, and in the 
collision and interaction of the vortex pairs (see the sec- 
ond and third pairs of panels in Fig. [5]), the vortices keep 
their identities and change partners, resulting into two 
new vortex pairs which propagate obliquely to the origi- 
nal propagation direction. On the other hand, as shown 

> 1) 



vortices that are multiply charged ( 
)le. Here the system of Eqs. (57) and (|58|) was 



in Fig. [6 

are unsta 

again solved numerically with the same initial condition 
as the one in Fig. [5j but with doubly charged vortices 
,si = +2, s-2 = —2, s 3 = —2, and S4 = +2. The second 
row of panels in Fig. [6] reveals that the vortex pairs keep 
their identities for some time, while a quasi ID density 
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FIG. 6 The electron density |*| 2 (left 
plot of the electron current i (\&V\E'* — 
associated with double charged (|s| = 2) 
obtained from a simulation of the time- 
Eqs. J571 and J581, at times t = 0, t 
t = 9.9 (upper to lower panels), with 



charged vortices dissolve into nonlinear 
turbulence. After |Shukla and Eliasson| ([ 



panel) and an arrow 
**V*) (right panel) 
2D electron vortices, 
dependent system of 
= 3.3, t — 6.6 and 
A = 5. The doubly 
structures and wave 
2006] ). 



cavity is formed between the two vortex pairs. At a later 
stage, the four vortices dissolve into complicated nonlin- 
ear structures and wave turbulence. Hence, the nonlinear 
dynamics is very different between singly and multiply 
charged solitons, where only singly charged vortices are 
long-lived and keep their identities. 



VI. QUANTUM ELECTRON FLUID TURBULENCE 

The statistical properties of quantum electron fluid 
turbulence and its associated electron transport proper- 
ties at nanoscales in a quantum plasma have been inves- 
tigated in both 2D an d 3D by using the coupled NL S 
and Poisson equations (Shaikh and Shukla 2007 2008 1. 



It has been found that nonlinear couplings between the 
EPOs of different scale sizes give rise to small-scale elec- 
tron density structures, while the ES potential cascades 
towards large-scales. The total energy associated with 



the quantum electron plasma wave turbulence processes 
a non-universal spectrum that depends on the quantum 
electron coupling parameter. 

To investigate 3D quantum electron plasma wave tur- 
bulence, we use the NLS- Poisson equations (Manfredi 
and Haas[ [200Tj |Shaikh and Shukla[ [2008l |Shukla[ |2006[ 



Shukla and Eliasson 2006) 



#V 2 W + - |*| 4/3 * = 0, 



and 



l*| 2 



1. 



(69) 



(70) 



were used, which govern the dynamics of nonlinearly in- 
teracting EPOs of different wavelengths. In Eqs. (69) 
and ( 70 1 the wavefunction is normalized by tJuq , the 



ES potential by fcsTp e /e, the time t by the electron 
plasma period w~ , and the space r by the electron 
Fermi-Thomas screening length A^e = VFe/^pe- Here 
y/~H = Tiujpe/ \i2kBTp e was introduced. 

The nonlinear wave-wave coupling studies have been 
performed to investigate the multi-scale evolution of a 
decaying 3D electron plasma wave turbulence, which is 
described by Eqs. (69) and (70). All fluctuations are ini- 



tialized isotropically (no mean fields are assumed) with 
random phases and amplitudes in Fourier space, and are 



evolved in time by the integration of Eqs. (69) and (70) 



numerically. The initial isotropic turbulent spectrum was 
initially chosen close to fc~ 2 , with random phases in all 
three directions. The choice of such (or even a flatter 
than k~ 2 ) spectrum treats the turbulent fluctuations on 
an equal footing and avoids any influence on the dynami- 
cal evolution that may be due to the initial spectral non- 
symmetry. 

The properties of 3D electron plasma wave turbu- 
lence, composed of nonlinearly interacting EPOs, were 
studied for two specific physical systems, corresponding 
to dense plasmas in the next generation of laser-based 



plasm a compression (LBPC) schemes ( Malkin et al. 
2007 ) , and in superdense astroph ysical objects (|Chabrier 



eTffi|20lffi|200l)i|Chabri^ and Lai[|2006 



Lai[ 2001 ) (e.g. white dwarfs). It is expected that in 

LBPC schemes, the electron number density may reach 

10 27 cm -3 and beyond. Hence, we have w pe = 1.76 x 10 18 

s" 1 , T Fe = 1-7 x 10~ 9 erg, hu} pe = 1.7 x 10~ 9 erg, and 

H = 1, and the electron Fermi-Thomas screening length 

\p e = 0.1 A. On the other hand, in the core of white 

dwarf stars, we typically have no ~ 10 30 cm~ 3 , yielding 
-1 t_ — 1 7 v 1 n-7 



u} pe = 5.64 x 10 19 s"\ T Fe = 1.7 x 10~ 7 erg (0.1 MeV), 
hu} pe = 5.64 x 10 



erg, H w 0.3, and X Fe = 0.025 A. The 
numerical solutions of Eqs. (69) and (70) for H = 0.4 
and H = 0.01 (corresponding to no = 10 27 cm~ 3 and 
n Q = 10 30 cm~ 3 , respectively) are displayed in Fig. [TJ 
which shows the electron number density and ES poten- 
tial distributions in the (x, y, z)-cube. 

Figure [7] reveals that the electron density distribution 
has a tendency to generate smaller length-scale struc- 
tures, while the ES potential cascades towards larger 
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FIG. 7 Small-scale fluctuations in the electron density re- 
sulting from steady turbulence simulations, for H = 0.4 (top 
panels) and H = 0.01 (bottom panels). Forward cascades are 
responsible for the generation of small-scale fluctuations seen 
in panels (a) and (c). Large scale structures are present in the 
ES potential, seen in panels (b) and (d), essentially resulting 
from an inverse cascade. After [Shaikh and Shukla] ( |2008[ ). 




FIG. 8 Energy per vector wavenumber of 3D EPOs in the 
forward cascade regime. A Kolmogorov-like spectrum Sk ~ 
fc -11 / 3 is observed for H = 0.4. The spectra l index changes 
as a function of H. After Shaikh and Shukla (20081. 



scales. The co-existence of the small and larger scale 
structures in turbulence is a ubiquitous feature of vari- 
ous 3D turbulence systems. For example, in 3D hydro- 
dynamic turbulence, the incompressible fluid admits two 
invariants, namely the energy and the mean squared vor- 
ticity. The two invariants, under the action of an external 
forcing, cascade simultaneously in turbulence, thereby 
leading to a dual cascade phenomena. In these pro- 
cesses, the energy cascades towards longer length-scales, 
while the fluid vorticity transfers spectral power towards 
shorter length-scales. Usually, a dual cascade is observed 
in a driven turbulence simulation, in which certain modes 



are excited externally through random turbulent forces 
in spectral space. The randomly excited Fourier modes 
transfer the spectral energy by conserving the constants 
of motion in fc-space. On the other hand, in freely decay- 
ing turbulence, the energy contained in the large-scale 
eddies is transferred to the smaller scales, leading to a 
statistically stationary inertial regime associated with the 
forward cascades of one of the invariants. Decaying tur- 
bulence often leads to the formation of coherent struc- 
tures as turbulence relaxes, thus making the nonlinear 
interactions rather inefficient when they are saturated. 
The power spectrum exhibits an interesting feature in 
the 3D electron plasma system discussed here, unlike the 



3D hydrodynam i c turb ulence ( Frisch 1995 |Kolmogorov 
1941a|bt |Lesiur[ p90| ). Fi gure [8] shows the energy £, 



per vector wavenumber. For isotropic 3D turbulence, it 
is related to the ener gy E/~ per scalar wavenumb er as 
[Knight and Sirovich| (l990[ ). For 



4?rfe £k\ see e.g. 



E, 

H = 0.4, the spectrum per vector wavenumber is close to 
£k ~ A;" 11 / 3 and hence yields the standa rd Kolmogorov 
power spectrum (Kolmogorov |l941a|b| ) E k ~ fc- 5 / 3 . 
However, the spectrum is not universal but changes for 
different values of H. For 2D quantum electron fluid 



turbulence ( Shaikh and Shukla 2007 1 the spectral slope 
was more c l ose to the Iroshn i kov-K raichnan power law 

k -3/2_ 



(Iroshnikov 1963[ |Kraichnan 
origin of the differences in the 



19G5) E, 



The 

observed spectral indices 
resides with the nonlinear character of the underlying 
plasma models, as nonlinear interactions in the 2D and 
3D systems are governed typically by different nonlinear 
forces. The latter modify the spectral evolution of tur- 
bulent cascades to a significant degree. Physically, the 
flatness (or deviation from the fc -5 / 3 law), results from 
the short wavelength part of the EPOs spectrum which 
is controlled by the quantum electron tunneling effect 
associated with the Bohm potential. The peak in the en- 
ergy spectrum can be attributed to the higher turbulent 
power residing in the EPO potential, which eventually 
leads to the generation of larger scale structures, as the 
total energy encompasses both the electrostatic poten- 
tial and electron density components. In the dual cas- 
cade process, there is a delicate competition between the 
EPO dispersions caused by the statistical pressure law 
(giving the k 2 Vp e term, which dominates at longer scales) 



and the quantum Bohm force (giving the h k /Am 2 term, 
which dominates at shorter scales). 

The electron diffusion in the presence of small and 
large scale turbulent EPOs can be estimated in the fol- 
lowing manner. An effective electron diffusion coefficient 
caused by the momentum transfer can be calculated from 
D eff = J °°(P(r,t) • P(r,t + t'))dt', where P is elec- 
tron momentum and the angular bracket denotes spa- 
tial averages and the ensemble averages are normalized 
to unit mass. The effective electron diffusion coefficient, 
D e ff, essentially relates the diffusion processes associ- 
ated with random translational motions of the electrons 
in nonlinear plasmonic fields. To measure the turbulent 
electron transport that is associated with the turbulent 
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FIG. 9 Time evolution of the effective electron diffusion coef- 
ficient associated with the large-scale ES potential and the 
small-scale electron density, for H = 0.4, H = 0.1 and 
H — 0.01. Smaller values of H corresponds to a small effec- 
tive diffusion coefficient, which characterizes the presence of 
small-scale turbulent eddies that suppress the electron trans- 
port. After IShaikh and Shuklal (12008 1. 



structures, -D e // is computed. It is observed that the 
effective electron diffusion is lower when the field pertur- 
bations are Gaussian. On the other hand, the electron 
diffusion increases rapidly with the eventual formation of 
large-scale structures, as shown in Fig. [9j The electron 
diffusion due to large scale potential distributions in a 
quantum plasma dominates substantially, as depicted by 
the solid-curve in Fig. [9] Furthermore, in the steady- 
state, nonlinearly coupled EPOs form stationary struc- 
tures, and D e ff eventually saturates. Thus, remarkably 
an enhanced electron diffusion results primarily due to 
the emergence of large-scale potential structures. 



VII. NONLINEARLY COUPLED EM AND ES WAVES 



A. Stimulated Scattering Instabilities 

First, we present the governing equations for the HF- 
EM waves and the EM wave driven modified EPOs and 



pe 



IPOs. We have flStenflo and Shukla[ |2009[ ) 

d 2 

— - c 2 V 2 

for the HF-EM wave, 
d 2 „ 3 



2 n l a 

W pe— A 



0, (71) 



dt 2 



v 2 v 2 



Am 2 



n 



2m 2 c 2 



rV|A| : 



for the HF-EM wave pressure driven EPOs, and 



(72) 



\dt 2 



r 2 



V 2 + 



4m P nii 



ni 
"o 



2m P m,c 2 



VIAI 



(73) 

for the EM wave pressure driven modified IPOs without 
the ion thermal, ion viscoelastic relaxation, and ion vis- 
cosity effects. Here Ctf — (kBTFe/nii) 1 ^ 2 and ni(<SC n ) 
is a small perturbation in the electron number density. 
Following the standard procedure of the parametric 



instabilities ( 


Murtaza and Shukla 




1984 


Sharma and 


Shukla 1983| Shukla 2006 Shukla and Eliasson 2006 


Shukla et al. 


1981; 


Yu et al. 




1974 


), we can Fourier an- 



obtain the nonlinear dispersion relations 



2 -< 2 k 2 \A \ 2 



2m 2 c 2 



and 



e 2 <fc 2 |A | 
2to p to,c 2 



1 



1 



i 



1 



(74) 



(75) 



for the driven EPOs and IPOs, respectively, which admit 
SRS, SBS, and modulational instabilities of the HF-EM 
pump (with the amplitude A ) in a quantum plasma. 
Here, D± = ±2cj (w-k- V g )-k 2 c 2 , where V g = kc 2 /2w 
is the group velocity of the HF-EM pump wave with the 
frequency uj = (fc 2 c 2 -I- w 2 ) : ^ 2 , and 



ni = uj 2 

R ^pe 



k 2 V 2 e 



n 2 k 4 

Am 2 ' 



(76) 



We turn our attention to nonlinear interactions be- 
tween large amplitude EM and ES waves in a quantum 
plasma. 



Shukla and Stenflo (2006) considered 
ear couplings between large amplitude EM 



nonlin- 
and 



waves 

finite amplitude electron and ion plasma waves, and pre- 
sented nonlinear dispersion relations that exhibit stimu- 
lated Raman scattering (SRS), stimulated Brillouin scat- 
tering (SBS), and modulational instabilities. The work 



of Shukla and Stenflo ( 2006 ) has been further generalized 



by including thermal corrections to the ES waves ( Sten- 



flo and Shukla 2009 ) and relativistic electron mass vari- 



ations (Shukla and Eliasson 2007) caused by EM waves 



in an unmagnetized quantum plasma. 



and 



k 2 C? 



TF 



ft 2 fc 4 



(77) 



The growth rates of SRS and SBS instabilities (Shukla 
and Stenflo 2006) are, respectively, 



1r 



uj pe eK \A Q \ 



and 



IB 



2v / 2tJoTTRm e c' 

Lo pe eK\A \ 
2\f2ujoQ,Brn e miC 



(78) 



(79) 
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The present results of SRS and SBS instabilities will help 
to identify the electrostatic spectral lines that are en- 
hanced by the large amplitude HF-EM pump wave in a 
quantum plasma. 



B. Nonlinearly Coupled Intense EM and EPOs 

Let us now consider nonlinear interactions between 
an arbitrary large amplitude circularly polarized electro- 
magnetic (CPEM) wave and nonlinear EPOs that are 
driven by the relativistic ponderomotive force (Shukla 



and Yu 1984 Shukla et al. 1985) of the CPEM waves 



Such an interaction gives rise to an envelope of the CPEM 
vector potential Aj_ = A±(x + iy) exp(— iw Q t + ikpz), 



which obeys the NLS equation ( jShukla and Eliasson 
20071) 



2ie 



d_ 

Of 



U, 



d_ 

dz 



d 2 A A 



1 



A ± + ~ ( — - 1 J A± = 0. 

_ (80) 

where e = uJo/uj pe , and the normalized (by */no) electron 
wavefunction ijj and the normalized (by moc 2 /e) scalar 
potential are governed by, respectively, 



in a standard manner ( Shukla et al.\ 1985) by letting 
<t>(z,t) = <f>t(z,t), A±(z,t) = [A + Ai(z,t)\exp(-ia t) 
and ip(z,t) = [1 + ip\(z, t)] exp(— i/3 i), where A is 
the large-amplitude CPEM pump and A\ is the small- 
amplitude perturbation of the CPEM wave amplitude 
due to the nonlinear coupling between the CPEM waves 
and EPOs, i.e. \A X \ < \Aq\, and ip x (< 1) is the 
small-amplitude perturbation in the electron wave func- 
tion. The constants ao and /3q are constant frequency 
shifts, determined from Eqs. (80) and (81) to be 
"o = (l/7o ~ l)/(2e), and (5 = (1 - -f )/H e , where 
70 = (1 + l^lol 2 ) 1 ^ 2 - The first-order perturbations in the 
electromagnetic vector potential and the electron wave 
function are expanded into their respective sidebands as 
A x {z,t) = A + exp(iKz - iftt) + A_exp(-iKz + iftt) 
and V'i( z j t) = i } + Gxp(iKz — iflt) + ip- exp(—iKz + iflt), 
while the potential is expanded as 4>{z, t) — <p exp{iK z — 
iflt) + <fi* exp(—iKz + iflt), where SI and K are the nor- 
malized frequency and the normalized wave number of 
the EPOs, respectively. Inserting the above mentioned 
Fourier ansatz into Eqs. ( 80 1— ( 82 ), linearizing the resul- 



tant system of equations, and sorting into equations for 
different Fourier modes, one obtains the nonlinear dis- 



persion relation ( Shukla and Eliasson , 2007 ) 



iH, 



chp H 2 d 2 j) 



+ (^-7 + 1)^ = 0, (81) 



and 




Dl) 2 7o 3 ' 



(83) 



d 2 4> 
dz 2 



|^| 2 -1, (82) 



where D ± = ±2e(n - KU g ) - K 2 and D L = 1 - e 2 + 
HlK 4 /4. One notes that Dl = yields the linear dis- 
persion relation f2 2 = 1 + H~l K 4 /4: for the EPOs in a 
H e = huj pe /moc is the ration between the plasmonic dense quantum plasma (JpTf^l pggTl. For H e 



where mo is the rest mass of the electrons, U g — koc/2ojQ 



energy density to the rest electron energy, and 7 
(1 + _l I 2 ) 1 / 2 is the relativistic gamma factor due to 
the electron quiver velocity in the CPEM wave fields. 
The time and space variables are in units of the elec- 
tron plasma period (<jJ~^) and the electron skin depth 
A e = c/u! pe . The electron density and A± a re in units 
of no and moc 2 /e Shukla and Eliassonj (2007). The non- 
linear coupling between intense CPEM waves and EPOs 
comes about due to the nonlinear current densit y, w hich 
is represented by the term |i/>| 2 Aj_/7 in Eq. (80). In 
Eq. 
tial 



we 

recover from (jS3h the nonlinear dispersion relation for 



relativistically large amplitude EM waves in a classical 
electron plasma (McKinstrie and Bingham 1992 1. The 



dispersion relation ( 83 1 governs stimulated Raman back- 



8li 1 



7 is the relativistic ponderomotive poten- 

The 



Shukla and Yu 1984 Shukla et al. 1985 



latter arises from the averaging (over the CPEM wave 
period 2it/luq) of the relativistic advection and the non- 
linear Lorentz force involving the electron quiver velocity 
and the CPEM wave electric and magnetic fields. 

A relativistically strong EM wave in a classical 
electron-ion plasma is subject to SRS and modulational 
instabilities (McKinstrie and Bingham 1992). One can 



ward and forward scattering instabilities, as well as the 
modulational instability. In the long wavelength limit 
Ug <C 1, e « 1 one can use the ansatz SI = ir, where the 
nor maliz ed (by u> pe ) growth rate I< 1, and obtain from 
Eq. H the growth rate T = (l/2)| J ftT|{(|A | 2 /7^)[l + 
K 2 /(T+ H 2 K 4 /4)} - K 2 } 1 / 2 of the modulational insta- 
bility. For I if I < 1 and H e < 1, the linear growth rate is 
only weakly depending on the quantum parameter H e . 

The quantum dispersion effects on nonlinearly cou- 
pled CPEM and EPOs can be studied by considering 
a steady state structure moving with a constant speed 



U g . Inserting the ansatz Aj_ = W(£) exp(— ift e t), ip = 



P(£) exp(ik e x 
where £ = z 



iuj e t) and <fi = </>(£) into Eqs 



Ugt 



U g /H e and 




(82), 



and 



expect that these instabilities will be modified at quan- 
tum scale by the dispersion effects caused by the tunnel- 
ing of electrons through the quantum Bohm potential. 
The growth rate of the relativistic parametric instabili- 
ties in a dense quantum plasma in the presence of a rel- 
ativistically strong CPEM pump wave can be obtained 



where W(£) and P(£) are real, one obtains from (f80|) - (f82j) 



the coupled system of equations (Shukla and Eliasson 



2007) 



d 2 W 



X - 



P 2 



1 )w = 0, 



(84) 
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Hi d 2 P , x 
where 7 = (1 + W 2 ) 1/2 , and 



7 + 1)^ = 0, 



P 2 



(85) 



(86) 



with the boundary conditions W = $ = and P 2 = 1 at 



frequency shift of the CPEM wave. In the limit H e — > 0, 
one has from (85) = 7—1, where P / 0, and one re- 



covers the classical (non-quantum) case of the relativistic 



solitary waves in a cold plasma (Marburger and Tooper 
1972b. 



The system of equations (84)— (86) admits a Hamilto- 
nian 



>h = 



_ 1 fdw\ 2 ^l(^_\ 2 _\(^± 



2 V dt, J 



2 \d£ 



1 



(87) 



(A + l)W 2 + P 2 - 7P 2 + 0P 2 - 4> = 0, 



0, w = 4> = 



where the boundary conditions d/dS, 
and |P| = 1 at |£| = 00 have been used. 

Numerical solutions of the quasi-stationary system 
-(86) are presented in Fi gs. [T75] an d 11 while time- 
dependent solutions of Eqs. (|80[)-(82 ) are displayed in 
Figs. [I2] and [13] Here parameters were used that 
are representative of the next generation of laser-based 



plasma compression (LBPC) schemes (Azechi et al. 



2006 



1985 



Malkin et al. 



2007]. The formula ( |Shukla et al 
10 \ S V1 will determine the nor 



_ eA ± /mc 2 = 6x10- 
malized vector potential, provided that the CPEM wave- 
length A s (in microns) and the CPEM wave intensity / 
(in W/cm 2 ) are known. It is expected that in LBPC 
schemes, the electron number density no may reach 10 27 
cm -3 and beyond, and the peak values of eA±/mc 2 
may be in the range 1-2 (e.g. for focused EM pulses 
with A s ~ 0.15 nm and J ~ 5 x 10 27 W/cm 2 ). For 

s _1 , one has huj pe = 1.76 x 10 -9 erg 
-7 



and H e = 0.002, since mc 2 = 8.1 x 10 -7 erg. The elec- 
tron skin depth A e ~ 1.7 A. On the other hand, a higher 
value of H e = 0.007 is achieved for u> pe = 5.64 x 10 18 s _1 . 
Thus, the numerical solutions below, based on these two 
values of H e , have focused on scenarios that are relevant 
for the next generation intense laser-solid density plasma 
interaction experiments (Malkin et al. 20071. 

Figures [lO] and [TT] exhibit numerical solutions of Eqs. 
d84b— dS6b for H e = 0.002 and H e = 0.007. The nonlinear 



boundary value problem was solved with the boundary 
conditions W = <fi = and P = 1 at the boundaries at 
£ = ±10. The solitary envelope pulse is composed of a 
single maximum of the localized vector potential W and 
a local depletion of the electron density P 2 , and a lo- 
calized positive potential <f> at the center of the solitary 
pulse. The latter has a continuous spectrum in A, where 
larger values of negative A are associated with larger am- 
plitude solitary EM pulses. At the center of the solitary 



EM pulse, the electron density is partially depleted, 
in panels a) of Fig. 



as 



10 and for larger amplitudes of 



the EM waves one has a stronger depletion of the clcc- 



10 



tron density, as shown in panels b) and c) of Fig 
For cases where the electron density goes to almost zero 
in the classical case (Marburger and Tooper, 1972), one 



|£| = 00. In Eq. (84), A = 2efi e represents a nonlinear 



important quantum effect is that the electrons can tun- 
nel through the depleted density region. This is seen in 
Fig. 11 where the electron density remains nonzero for 
0.007 in panels a), while the density shrinks to zero 



for H e = 0.002 in panel b). 



a) X=-0.3 



b) ^=-0.34 
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FIG. 10 The profiles of the CPEM vector potential W (top 
row), the electron number density P 2 (middle row) and the 
scalar potential $ (bottom row) for A = —0.3 (left column), 
A = —0.34 (middle column) and A = —0.4 (right column), 
with H e = 0.002. After iShukla and Eliasson|(|2007|. 



Figures 12 and 13 depict numerical simulation results 
of Eqs. ( 80 )-( 82 ) for the long-wavelength limit character- 



ized by Wo ~ 1 and V g ~ 0. As initial conditions, we used 
an EM pump with a constant amplitude Aj_ = Aq = 1 
and a uniform plasma density ip = 1, together with a 
small amplitude noise (random numbers) of order 10 -2 
added to A± to give a seeding any instability. The 
numerical results are displayed in Figs. [12] and [13] for 
H e = 0.002 and H e = 0.007, respectively. In both cases, 
we see an initial linear growth phase and a wave collapse 
at t a 70, in which almost all the CPEM wave energy 
is contracted into a few well separated large-amplitude, 
localized CPEM envelopes, associated with an almost 
complete depletion of the electron density at the center 
of the CPEM wavepacket, and a large- amplitude posi- 
tive electrostatic potential. One can see that there is a 
more complex dynamics of localized CPEM wavepackets 
for H e = 0.007, shown in Fig. [13] in comparison with 
H e = 0.002, shown in Fig. [12] where the wavepackets 
are almost stationary when they are fully developed. 
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a) He=0.007 



b) He=0.002 







-5 5 -5 5 

FIG. 11 The profiles of the CPEM vector potential W (top 
row), the electron number density P 2 (middle row) and the 
scalar potential $ (bottom row) for H e = 0.007 (left column) 
and H e = 0.002 (rig ht column), with A = -0.34. After |Shuk la 
and Eliasson] ( |2007[ ). 
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FIG. 12 The dynamics of the CPEM vector potential A± 
and the electron number density \4>\ 2 (upper panels) and of 
the electrostatic potential $ (lower panel) for H c — 0.002. 
After IShukla and Eliassonl (|2007h . 



VIII. MAGNETIZED QUANTUM PLASMAS 

Magnetized quantum plasmas occur in white dwarf 
stars and on the surface of magnetized stars (e.g. 
magnetars) where degenerate electrons could be ultra- 
relativistic, but the ions are in a non-degenerate state. 
How strong magnetic fields in dense stars come about is 
still a mystery, although there are evidence of the strong 
magnetization of dense plasmas in astrophysical environ- 
ments. In dense magnetized plasmas, one has to account 
for the Lorentz force and the Landau quantization effect 
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FIG. 13 The dynamics of the CPEM vector potential A± 
and the electron number density \ip\ 2 (upper panels) and the 
electrostatic potential <j> (lower panel) for H e — 0.007. After 
Shukla and Eliassonl 12007b. 



(Landau and Lifshitz, 1998a), and develop the appropri 



ate quantum Hall-magnetohydrodynamics (Q-HMHD) 
equations starting from the Wigner-Maxwell equations. 
We stress, however, that the Q-HMHD equations dis- 
cussed here does not capture the particular physics of 
the quantized Hall resistance Rk — fi/ve 2 (Klitzing et 
al. 1980). In semiconductors with 2D electrons, the 



latter is associated with the quantized electron density 
n q = vcBq/Tic at high magnetic fields and low temper- 
ature, where v is an integer, appearing in the electron 
current (— en q Ud) flowing through a conductor. Here 
Ud = (c/£>q)E x B is the cross-field electron drift asso- 
ciated with the space charge electric field E that results 
from the motion of electrons by the Lorentz force. The 
Ohm's law, in turn, determines the von Klitzing resis- 
tance, which is independent of the magnetic field. 



A. Landau Quantization 

In a strong magnetic field zBq, where z is the unit 
vector along the z-axis in a Cartesian coordinate system, 
and Bo the strength of the external magnetic field, the 
electron motion in a plane perpendicular to the mag- 



netic field direction is quantized (Landau and Lifshitz 



1998b). The electron energy level is determined by the 



non-relativistic limit by the expression (Landau and Lif- 



shitz 1998b Tsintsadze 2010) 



e. 



La 



2m e 



(21 + 1 + a)n B B , 



where p z is the electron momentum in the z— direction, 
I the orbital angular number (I = 0,1,2), and a = ±1 
represents the spin orientation. For a = —1, we have 
from 
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2m e 



lhuj r 



(89) 



where uj ce = eBo/m e c is the electron gyrofrequency. 
Accordingly, the Fermi-Dirac electron distribution is 
dTsintsadzel[2010) 



F D (p z ,l) oc 



1 + CXp [(E z + lTlW ce - fi e ) /k B T e 



(90) 



where E z = (m e /2)v 2 is the parallel (to z) kinetic energy 
of degenerate electrons. 

Assuming that \lio ce — (i e \ ^ k B T e , one can ap- 
proximate the Fermi-Dirac distribution function by the 
Heaviside step function H(p e — £y), which equal 1 for 
p, e = E Fe = k B T Fe = (p F /2m e ) 1 / 2 > £ l e and zero for 
Ep e < £ e I where p F = m e VTF- The equilibrium electron 
number density is ( Tsintsadze| |201Q ) 



P 3 f 
2tt 2 Ti 3 



n.3/2 



(91) 



where Tg = fiuj ce /k B T Fe . The current carried by 
degenerate electrons in a magnetized quantum plasma 
is —en e Ud, which yields the plasma resistivity R s — 
en e c/B Q . 




FIG. 14 Dispersion curves for EB waves in a Fermi-Dirac 
distributed plasma, showing several EB modes and the UH 
cutoff. After |Eliasson and Shukla| ( |2008b| ) . 



that the propagating UH waves have positive (negative) 
group dispersion in plasmas with ui pe > (<)v3w ce . 

Furthermore, the refractive index N x for the EP-EM 
waves propagating along the x axis (which is orthogonal 
to z) is (|Shukla| [2007|> 



B. ESOs and EM Waves 

In a magnetized quantum plasma, there are finite den- 
sity perturbations associated with high-frequency elec- 
trostatic electron-Bernstein (EB) waves and elliptically 
polarized EM waves (EP-EM waves) that propagate 
across the magnetic field direction z. Furthermore, the 
CPEM wave propagating along z are not associated with 
any density perturbation. 

The dispersion relation for the EB waves in a Fermi- 
Dirac distributed plasma is in the ultra-cold limit ( Elias 
son and Shuklal 12008b I 



3^ 2 



dO 



sin(f20) sin(0) sin(£ e ) — £ e cos(£ e 



= 0, 
(92) 



where ft = uj/uj ce , £ e = (2k\p Fe ) cos(0/2), and p Fe = 
VF e /ujp e is the gyroradius of degenerate electrons. So- 
lutions of Eq. (92) are plotted in Fig. [Ml for the case 
^UH = 4cj ce , where Uuh — (^p e + "L) 1 "' 2 is the upper- 
hybrid (UH) resonance frequency. In the long wavelength 
limit (viz. k\p Fe < 1), Eq. (92) yields 



3 ^lAyh 



5 (w 



(93) 



where k± is the perpendicular (to z component of the 
propagation wave vector. For u sa ojh, Eq. (93) reveals 



N,. 



= 1 - 



, , 2 i , 2 



l + V (a)k 2 x \l] 



(94) 



J UH 



where k x is the x component of the propagation wave 
vector, A 2 = h 2 /4m e V Fe , X b = y / h/2m e ui ce , rj(a) = 
2tanh(a), a = p, B Bo/k B T Fe . Several comments are 
in order. First, we note that the electron spin-1/2 ef- 
fect enhances the electron gyrofrequency by a factor 
(1 + rik'l.X 2 ) 1 / 2 in the numerator of the third term in the 
right-hand side of (94). Second, the quantum Bohm force 
produces a dispersion term Tik A /Am 2 in the denominator 
of the third term in (94). Third, in the limit of vanishing 
Ti, Eq. (95) correctly reproduces the EP-EM wave dis- 
persion relation. Furthermore, Eq. (94) reveals that the 
cut-off frequencies (at k x = 0) in dense magnetoplasmas 
are 



{4u, 



2 

pe 



,2 



(95) 



which are the same as the cutoffs of the X (upper sign) 
and Z (lower sign) mode waves in a classical plasma 
(Chen) 2006). Short wavelength electromagnetic propa- 



gation in magnetized quantum plasmas, including quan- 
tum electrodynamic effects, has also been considered by 



Lundin et al. (2007) 



The vector representation of spinning quantum par- 
ticles in the quantum theory was first introduced by 
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Takabayasi ( 1955 ) who developed the QHD involving 
the evolution of the quantum particle spin. The idea 
of Takabayasi has been further elaborated by Brodin et 
al. (2010) in the context of the spin contribution to the 



ponderomotive force of the magnetic field-aligned CPEM 
waves in a quantum magnetoplasma. In fact, by using the 
non-relativistic electron momentum equation (Brodin et 



al. 2010) 



in, [ — + u e • V ) u e 



1 



e E + u e x B 



|mbV(B-s), 



(96) 



and the spin evolution equation 

d „\ gpB 



dt 



u e • V s 



(B x s) . 



(97) 



together with Ampere's law and suitable Maxwell's equa- 
tion (incorporating the electron magnetization current, 
Jm = — (47r/c)(g^s/?i)V x (n e x s), due to the electron 
1/2— spin effect), where s is the spin angular momentum, 
with its absolute value |s| = sq = fi/2. The quantity 
g = 2.0023192 is the electron Gaunt factor (sometimes 
called the g factor or spectroscopic splitting factor) . The 
value g = 2 is predicted from Dirac's relativistic theory of 
the electron, while the correction to this value comes from 



the quantum electrodynamics (Bransden and Joachain 
~ T996I 



2000 Kittel 



Brodin et al. (20101 derived the spin ponderomotive 
force zF for the CPEM wave, where 



F 



2 ? 

. STUb 



sq 



mlh (w ± u) g ) 



d 



d 



dz (u) ± ujg ) dt 



IB,, 



(98) 



Here u> g = gpsBo/ti the spin-precession frequency, and 
B^, is the CPEM wave magnetic field. The spin pondero- 
motive force comes from the averaging of the third term 
in (99) over the CPEM wave period 2-k/cj. The CPEM 
wave frequency uj is determined from the dispersion re- 
lation 



it 



(U) ± Wg) 



Ni 



uj(w ± w ce ) ' 



(99) 



where N z — k z c/uj, k z is the component of the wave vec- 
tor k along the z axis, uj^ = g 2 So/4m e Xl, A e = c/w pe , 
and the +(— ) represents the left- (right-) hand circular 
polarization. The w^-term in (100) is associated with 
the electron spin evolution. It changes the dispersion 
properties of the magnetic field-aligned EM electron- 
cyclotron waves in a quantum magnetoplasma. Fur- 
thermore, the spin-ponderomotive force induces a strong 
spin-polarization of a quantum magnetoplasma. 

It should be noted that there is also a standard non- 



stationary ponderomotive force (zF e ) (Karpman and 



Washimi 19771 of the CPEM waves arising from the av- 



eraging of the nonlinear Lorentz force term — (e/m e c)z ■ 
(u e x Bu,) over the CPEM wave period 2ir/ui, where 



F„ 



2m1ui{ui ± oj ce ) 



d 



± 



d 



dz ± uj ce ) dt 



|Et„| 2 , 



(100) 

where E„ = (uj/k z c)'B w is the CPEM wave electric field. 



C. Q-HMHD Equations 

To a first approximation, the dynamics of low phase 
speed (in comparison with the speed of light in vacuum) 
electromagnetic waves in dense magnetoplasmas is mod- 
eled by the Q-HMHD equations. The latter include the 
inertialess electron momentum equation 



E 



u P x B 



VF C , 



(101) 



where the quantum Bohm and quantum spin forces are 
supposed to be unimportant on the characteristic scale- 
length of present interest. The degenerate electrons are 
coupled with the non-degenerate ions through the EM 
forces. The ion dynamics is governed by the ion continu- 
ity equation (43) and the momentum equation 



mini 



dt 



E- 



1 



u, x B 



(102) 



where d/dt — (d/dt) +Uj • V. For the sake of simplicity, 
we have here assumed that T m d/dt <C 1 and dui/dt ^> 
{i)/pi)V ■ Vu; + pr 1 ^ + r//3)V(V • Ui). The EM fields 
are given by Ampere's law 



dB 

~dt : 

and Maxwell's equation 



-cV x E, 



^ ^ Aire, , IdE 

V x B = (rajUj - n e Ui) H — . 

c c dt 



(103) 



(104) 



By using (102), we can eliminate the electric field E 
from (103), obtaining for a quasi-neutral (n e = = n) 
quantum magnetoplasma 



rriin 



dui 
dt 



-VPc 



VB 2 



(B- V)B 



(105) 



where we have used (104) without the displacement cur- 
rent (the last term on the right-hand side) for the low- 
phase speed (in comparison with c) EM wave phenomena. 
By using the electric field from (102), we can write (103) 
as 
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SB 

~dt 



V x (uj x B) 



e dt 



(106) 



Equations (43), (105) and (106) are the desired Q- 
HMHD equations for studying the linear and nonlin- 
ear dispersive EM waves, as well as new aspects of 3D 
quantum fluid turbulence in a quantum magnetoplasma 
with degenerate electrons having Chandrasekhar's pres- 
sure law. However, when the Landau quantization effect 
in a very strong magnetic field is accounted fo r, one can 
rep lace P c by the appropriate pressure law (Eliezer et 
aTl[2005| 



Pt 



AeB (2m e ) 1 / 2 E 3 F /2 
3(27r) 2 n 2 c 

l + 2 £(l- 

i=i v 



k B T Fe 



3/2' 



(107) 



where the value of l m is fixed by the largest integer that 
satisfies ksTpe — lhcu ce < 0. 



IX. SUMMARY AND OUTLOOK 

In this Colloquium paper, we have described the essen- 
tial physics of quantum plasmas with degenerate electron 
fluids. We have reviewed the properties of quantum plas- 
mas and quantum models that describe the salient fea- 
tures of linear and nonlinear ES and EM waves. Specifi- 
cally, the focus of the present colloquium article has been 
on developing the model nonlinear equations that depict 
new features of nonlinear waves and quantum electron 
fluid turbulence at nanoscales. Numerical simulations of 
the nonlinear Schrddinger (NLS)-Poisson equations re- 
veal quasi-stationary, localized structures in the form 
of one-dimensional electron density holes (dark solitons) 
and 2D quantum electron vortices. These localized quan- 
tum structures, which are associated with a local deple- 
tion of the electron density and a positive electrostatic 
potential, arise due to a balance between the nonlinear 
and dispersion effects involved in the dynamics of nonlin- 
early interacting EPOs. In 2D, there exist a class of quan- 
tum electron vortices of different excited states (charge 
states). Furthermore, numerical simulations also depict 
that the time-dependent NLS-Poisson equations exhibit 
stability of a dark soliton in one-space dimension. In 2D 
the dark solitons of the first excited state are stable and 
the preferred nonlinear state is in the form of quantum 
vortex pairs of different polarities. The one-dimensional 
dark soliton and singly charged 2D quantum vortices are 
thus long-lived nonlinear structures at nanoscales. Also 
presented are theoretical and computer simulation stud- 
ies of nonlinearly coupled intense EM waves and EPOs 
in an unmagnetized quantum plasma. We have reported 
new classes of stimulated scattering instabilities of EM 



waves and trapping of intense EM waves in a quantum 
electron density hole. 

It should be noted that inclusion of non-degenerate ion 
dynamics gives rise to new features to linear and nonlin- 
ear IPOs (lEliasson and Shuklal |2008a| IHaas et aU [2003 1. 



Furthermore, nonlinear equations governing the coupling 
between the dispersive Langmuir and ion-acoustic waves, 
which are known as the quantum Zakharov equations 



(Garcia et al. 2005 Haas and Shukla 2009 Misra et 



al. 2008 Simpson et al. 20091, admit periodic, quasi- 



periodic, chaotic and hyper-chaotic states (Misra et al. 



20081, in addition to arresting the Langmuir wave col 



lapse (Haas and Shukla 2009 Simpson et al. 2009) due 



to quantum dispersion effects. There may also emerge 
new aspects of nonlinear EPOs and IPOs when the par- 
ticle trapping ( |Jovanovic and Fedele 2007 ) in the strong 
wave potential is included. Here one has to obtain non- 
linear solutions of non-stationary Wigner-Poisson equa- 
tions, which might reveal a modified (by the electrostatic 
wave potential) Fermi-Dirac electron distribution func- 
tion. Furthermore, there is a scope for studying the col- 
lective nonlinear response of correlated Coulomb electron 
systems at finite temperatures by means of kinetic the- 
ory concepts (Domps et al. 1997| ) to incorporate colli- 
sions and Green's function methods originally developed 
by Baym and Kadanoff ( 1961 ). We note that the Baym- 
Kadanoff approach has been used by |Kwong and Bonitz| 
(2000) to investigate the dielectric properties (viz. in- 



verse dielectric function and dynamic structure factor) 
of linear EPOs in a correlated electron gas. Further- 
more, the ion-ion dynamic structure factor, which con- 
tains a wealth of information about ions including struc- 
ture and low-frequency collective modes in a dense quan- 
tum plasma, has been studied by Murillo (2010). 



The field of the nonlinear quantum plasma physics is 
vibrant, and its potential applications rest on our com- 
plete understanding of numerous collective processes in 
compact astrophysical objects, as well as in the next- 
generation of intense laser-solid density plasma experi- 
ments and in the plasma assisted nanotechnology (e.g. 
quantum free electron laser devices, quantum-diodes, 
metallic nanostructures, nanowires, nanotubes, etc.). 
However, nonlinear quantum models presented in this 
Colloquium paper have to be further improved and gen- 
eralized by including the effects of the electron exchange 
interactions, strong electron-electron correlations, equi- 
librium inhomogeneities of the magnetic field and the 
plasma density, as well as fully relativistic and Lan- 
dau quantization effects in a nonuniform quantum mag- 
netoplasma. We have also to understand the features 
of quantum oscillations of electrons and possible forma- 
tion of bound states of electrons in the presence of an 
external magnetic field. For this purpose, we have to 
calculate the interaction potential among highly corre- 
lated electrons and use molecular dynamic simulations 
to demonstrate attraction among electrons due to collec- 
tive wave-quantum particle interactions that give rise to 
Cooper's pairing of degenerate electrons. Cooper's pair- 
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ing of electrons could possibly provide a scenario of super- 
conducting behavior of a quantum plasma. Furthermore, 
2D system composed of electron clusters at finite tem- 
perature exhibits Wigner Coulomb crystallization (Eg- 



ger et al. 1999 Filinov et al. 2001). The latter has 



been investigated by means of Monte Carlo simulations 
based on a quantum Hamiltonian with parabolic confin- 
ing and Coulomb interaction potentials. In a nonuniform 
quantum magnetoplasma, we ha ve ES drift waves (|Ali| 
~eTal\ [20071 ISaleem et aT| [20081 [Shokri and Rukhadzef 
1999) which can drastically affect the cross-field elec- 
tron transport. For applications to plasma assisted 
nano-technology devices (e.g. nonlinear electrostatic and 
electromagnetic surface waves in metallic nanostructure- 
devices, photonic band gap and x-ray optical systems, 
quantum X-ray free-electron laser systems), one must 
also study nonlinear collective processes by including 
both the electron spin-1/2 and quantum electron tunnel- 
ing effects on an equal footing. Finally, the localization 
of coupled ES and EM waves due to nonlinear quantum 
effects in a nonuniform quantum magnetoplasma with 
an arbitrary electron pressure degeneracy should provide 
clues to the origin of v ery intense X-rays (|Coe et al.\ 
1978) and gamma rays (Hurley et al. 20051 from both 



astrophysical and laboratory plasmas. 



Acknowledgments 

This research was supported by the Deutsche 
Forschungsgemeinschaft through the project SH21/3-1 of 
the Research Unit 1048, and by the Swedish Research 
Council (VR). 



References 

Ali, S., N. Shukla, and P. K. Shukla, 2007, Europhys. J. Lett. 
78, 45001. 

Adolfath, R. M., A. G. Petukhov, and I. Zutic, 2008, Phys. 

Rev. Lett. 101, 207202. 
Anderson, D., B. Hall, M. Lisak, and M. Marklund, 2002, 

Phys. Rev. E 65, 046417. 
Andreev, A. V., 2000, JETP Lett. 72, 238. 
Ang, L. K., T. J. T. Kwan, and Y. Y. Lau, 2003, Phys. Rev. 

Lett. 91, 208303. 
Ang, L. K., and P. Zhang, 2007, Phys. Rev. Lett. 98, 164802. 
Ancona, M. G., and G. J. Iafrate, 1989, Phys. Rev. B 39, 

9536. 

Atwater, H. A., 2007, Sci. Am. 296, 56. 

Azechi, R, et al., 2006, Plasma Phys. Control. Fusion 48, 
B267. 

Balescu, R., and W. Zhang, 1988, J. Plasma Phys. 40, 215. 
Barnes, W., A. Dereux, and T. Ebbesen, 2003, Nature (Lon- 
don) 424, 824. 

Baym, G., and L. P. Kadanoff, 1961, Phys. Rev. 124, 287. 
Becker, K. H, K.H. Schoenbach, and J.G. Eden, 2006, J. 

Phys. D: Appl. Phys. 39, R55. 
Benvenuto O. G., and M. A. De Vito, 2005, Mon. Not. R. 

Astron. 362, 891. 



Berestetskii, B., E. M. Lifshitz, and L. P. Pitaevskii, 1999, 
Quantum Electrodynamics (Butterworth-Heinemann, Ox- 
ford), p. 123. 

Bohm, D., 1952, Phys. Rev. 85, 166. 

Bohm, D., 1953, Phys. Rev. 92, 626. 

Bohm, D., and D. Pines, 1953, Phys. Rev. 92, 609. 

Bonitz, M., 1998, Quantum Kinetic Theory (Teubner, 
Stuttgart). 

Bransden, B. H, and C. J. Joachain, 2000, Quantum Me- 
chanics (2nd Edition) (Pearson Education Limited, Essex, 
England). 

Brittin, W. E., and W. R. Chappell, 1962, Rev. Mod. Phys. 
34, 620. 

Brodin, G., and M. Marklund, 2007a, New J. Phys. 9, 277. 
Brodin, G., and M. Marklund, 2007b, Phys. Plasmas 14, 
11207. 

Brodin, G., and M. Marklund, 2007c, Phys. Rev. E 76, 
055403(R). 

Brodin, G., M. Marklund and G. Manfredi, 2008a, Phys. Rev. 

Lett. 100, 175001. 
Brodin, G., M. Marklund, J. Zamanian, A. Ericsson, and P. 

L. Mana, 2008b, Phys. Rev. Lett. 101, 245002. 
Brodin, G., A. P. Misra, and M. Marklund, 2010, Phys. Rev. 

Lett. 105, 105004. 
Burt, P., and D. Wahlquist, 1962, Phys. Rev. 125, 1785. 
Carruthers, P., and F. Zachariasen, 1983, Rev. Mod. Phys. 

55, 245. 

Chabrier, G., et al, 2002, J. Phys.: Condens. Matter 14, 
9133. 

Chabrier, G., D. Saumon, and A. Y. Potekhin, 2006, J. Phys. 
A: Math. Gen. 39, 4411. 

Chabrier, G., 2009, Plasma Phys. Control. Fusion 51, 124014. 

Chandrasekhar, S., 1931a, Astrophys. J. 74, 81. 

Chandrasekhar, S., 1931b, Phil. Mag. 11, 592. 

Chandrasekhar, S., 1935, Mon. Not. R. Astron. Soc. 170, 405. 

Chandrasekhar, S., 1939, An Introduction to the Study of Stel- 
lar Structure (University of Chicago Press, Chicago), p. 
360. 

Chang, D. E., A. S. S0rensen, P. R. Hemmer, and M. D. 
Lukin, 2006, Phys. Rev. Lett. 97, 053002. 

Chen, F. F., 2006, Introduction to Plasma Physics and Con- 
trolled Fusion. Volume 1, Plasma Physics., Second edition 
(Springer, New York). 

Coe, M. J., A. R. Engel, and J. J. Quenby, 1978, Nature 
(London) 272, 37. 

Cowley, S. C, R. M. Kulsrud, and E. Valeo, 1986, Phys. 
Fluids 29, 430. 

Crouseilles, N., P. A. Hervieux, and G. Manfredi, 2008, Phys. 
Rev. B 78, 155412. 

Dharma-wardana, C, and F. Perrot, 1995, in Density Func- 
tional Theory, Eds. E. K. U. Gross and R. M. Dreizler 
(Plenum Press, New York). 

Dirac, P. A. M., 1981, Principles of Quantum Mechanics (Ox- 
ford University Press, Oxford). 

Domps, A., P.-G. Reinhard, and E. Suraud, 1997, Ann. Phys. 
26, 171. 

Drake, R. P., 2009, Phys. Plasmas 16, 055501. 
Drake, R. P., 2010, Phys. Today 63, 28. 

Egger, R., W. Hausler, C. H. Mak, and H. Grabert, 1999, 

Phys. Rev. Lett. 82, 3320. 
Eliasson, B., and P. K. Shukla, 2010, Relativistic 

Laser-Plasma Interactions in the Quantum Regime, 

arXiv:1011.5801 vl [physics. plasm-ph]. 
Eliasson, B., and P. K. Shukla, 2008a, J. Plasma Phys. 64, 



23 



581. 

Eliasson, B., and P. K. Shukla, 2008b, Phys. Plasmas 15, 
102102. 

Eliezer, S., P. Norreys, J. T. Mendonga, and K. Lancaster, 

2005, Phys. Plasmas 12, 052115. 
Else, D., R. Kompaneets, and S. V. Vladimirov, 2010, Phys. 

Rev. E 82, 026410. 
Ferry, D. K., and J.-R. Zhou, 1993, Phys. Rev. B. 48, 7944. 
Feynman, R. P., and H. Kleinert, 1986, Phys. Rev. A 34, 

5080. 

Filinov, A. V., M. Bonitz, and Yu. E. Lozovik, 2001, Phys. 

Rev. Lett. 86, 3851. 
Fortov, V. E., 2009, Phys. Usp. 52, 615. 

Frisch, U., 1995, Turbulence (Cambridge University Press, 
Cambridge) . 

Froula, D. H., S. H. Glenzer, N. C. Luhmann, Jr., J. 
Sheffield, 2011, Plasma Scattering of Electromagnetic Ra- 
diation: Theory and Measurement Techniques, 2nd Edition 
(Academic Press, New York). 

Garcia, L. G., F. Haas, L. P. L. de Oliveira, and J. Goedert, 
2005, Phys. Plasmas 12, 012302. 

Gardner, C. L., and C. Ringhofer, 1996, Phys. Rev. E 53, 
157. 

Ghoshal, A., and Y. K. Ho, 2009a, Phys. Rev. A 79, 062514. 
Ghoshal, A., and Y. K. Ho, 2009b, J. Phys. B: At. Mol. Opt. 

Phys. 43, 175006. 
Glenzer, S. H., and R. Redmer, 2009, Rev. Mod. Phys 81, 

1625. 

Glenzer, S. H., et al, 2007, Phys. Rev. Lett. 98, 065002. 
Gottlieb, D., and S. A. Orszag, 1977, Numerical Analysis of 

Spectral Methods (SIAM, Philadelphia) . 
Gregori, G., and D. O. Gericke, 2009, Phys. Plasmas 16, 

056306. 

do Groot, S. R., and L. G. Suttorp, 1972, Foundations of 

Electrodynamics (North- Holland, Amsterdam). 
Guillot, T., 1999, Science 286, 72. 

Gursky, H., 1976, in Frontiers of Astrophysics, Ed. E. H. 

Avrett (Harvard University Press, Cambridge, Massachus- 

sets), Chap. 5, pp. 152,153. 
Haas, F., G. Manfredi, and M. R. Feix, 2000, Phys. Rev. E 

62, 2763. 

Haas, F., 2005, Phys. Plasmas 12, 062117. 

Haas, F., 2007, Europhys. Lett. 44, 45004. 

Haas, F., L. G. Garcia, J. Goedert, and G. Manfredi, 2003, 
Phys. Plasmas 10, 3858. 

Haas, F., and P. K. Shukla, 2009, Phys. Rev. E 79, 066402. 

Haas, F., M. Marklund, G. Brodin, and J. Zamanian, 2010a, 
Phys. Lett. A 374, 481. 

Haas, F., J. Zamanian, M. Marklund, and G. Brodin, New. 
J. Phys. 12, 073027 (2010b). 

Hakim, R., and J. Heyvaerts, 1978, Phys. Rev. A 18, 1250. 

Harding, A. K., and D. Lai, 2006, Rep. Prog. Phys. 69, 2631. 

Haug, H., and S. W. Koch, 2004, Quantum Theory of Op- 
tical and Electronic Properties of Semiconductors (World 
Scientific, Singapore). 

Haug, H., and A. -P. Jauho, 2007, Quantum Kinetics in Trans- 
port and Optics of Semiconductors, Springer Series in Solid- 
State Sciences, Vol. 123 (Springer, Berlin). 

Hillery, M., et al, 1984, Phys. Rep. 106, 121. 

Hohenberg, P., and W. Kohn, 1964, Phys. Rev. 136, B864. 

Holland, P. R., 1993, The Quantum Theory of Motion (Cam- 
bridge University Press, Cambridge). 

Horn, H. M., 1991, Science 252, 384. 

Hu, S. X., and C. H. Keitel, 1999, Phys. Rev. Lett. 83, 4709. 



Hurley, K., S. E. Boggs, D. M. Smith et al, 2005, Nature 
(London) 434. 

Iafrate, G. J., H. L. Grubin, and D. K. Ferry, 1981, J. Phys. 

(Paris) (Colloq.) 42, C10-307. 
Ichimaru, S., 1982, Rev. Mod. Phys. 54, 1017. 
Ichimaru, S., and S. Tanaka, 1986, Phys. Rev. Lett. 56, 2815. 
Iroshnikov, P., 1963, Sov. Astron. 7 566. 
Jovanovic, D., and R. Fedele, 2007, Phys. Lett. A 364, 304. 
Jiingel, A., D. Matthes, and J. F. Milisic, 2006, SIAM 67, 46. 
Hohenberg, P., and W. Kohn, 1964, Phys. Rev. 136, 864. 
Kaw, P. K., and A. Sen, 1998, Phys. Plasmas 5, 3552. 
Karpman, V. I., and H. Washimi, 1977, J. Plasma Phys. 18, 

173. 

Kelly, D. C, 1964, Phys. Rev. 134, A641. 

Kittel, C, 1996 Introduction to Solid State Physics, Seventh 

Edition (John Wiley & Sons, Inc., New York) 
Kleinert, H., 1986, Phys. Lett. B 181, 324. 
Klimontovich, Y. L., and V. P. Silin, 1952a, Doklady Akad. 

Nauk. SSSR 82, 361. 
Klimontovich, Y. L., and V. P. Silin, 1952b, Zh. Eksp. Teor. 

Fiz. 23, 151. 

Klimontovich, Y. L., and V. P. Silin, 1961, in Plasma Physics, 
Ed. J. E. Drummond (McGraw Hill, New York), pp. 35-87. 

v. Klitzing, K., G. Dorda, and M. Pepper, 1980, Phys. Rev. 
Lett. 45, 494. 

Knight, B., and L. Sirovich, 1990, Phys. Rev. Lett. 65, 1356. 
Koester, D., and G. Chanmugam, 1990, Rep. Prog. Phys. 53, 
837. 

Kohn W., and L. J. Sham, 1965, Phys. Rev. 140, A1133. 
Kolmogorov, A. N., 1941a, Dokl. Akad, Nauk SSR 30, 301. 
Kolmogorov, A. N., 1941b, Dokl. Akad, Nauk SSR 31, 438. 
Kraichnan, R. H., 1965, Phys. Fluids 8, 1385. 
Kremp, D., Th. Bornath, M. Bonitz, and M. Schlanges, 1999, 

Phys. Rev. E 60, 4725. 
Kremp, D., M. Schlanges, and W. D. Kraeft, 2005, Quantum 

Statistics of Nonideal Plasmas (Springer, Berlin). 
Kritcher, A. L., P. Neumayer, J. Castor et al, 2008, Science 

322, 69. 

Kuzelev, M. V., and A. A. Rukhadze, 1999, Phys. Usp.42, 
603. 

Kwong, N.-H., and M. Bonitz, 2000, Phys. Rev. Lett. 84, 
1768. 

Lai, D., 2001, Rev. Mod. Phys. 73, 629. 

Landau, L. D., andE. M. Lifshitz, 1998a, Quantum Mechanics 

(Butterworth-Heinemann, Oxford) . 
Landau, L. D., and E. M. Lifshitz, 1998b Statistical Physics 

(Butterworth-Heinemann, Oxford) . 
Lau, Y. Y., D. Chernin, D. G. Colombant, and P. T. Ho, 1991, 

Phys. Rev. Lett. 66, 1446. 
Lee, H. W., 1995, Phys. Rep. 259, 147. 

Lee, H. J., P. Neumayer, J. Castor, et al, 2009, Phys. Rev. 

Lett. 102, 115001. 
Lesieur, M., 1990, Turbulence in Fluids (Kluwer, Dordrecht). 
Lifshitz, E. M., and L. P. Pitaevskii, 1981, Physical Kinetics 

(Butterworth-Heinemann, Oxford), pp. 164-165. 
Lindhard, D. J., 1954, Mat. Fys. Medd. K. Dan. Vidensk. 

Selsk. 28, 8. 
Lindl, J., 1995, Phys. Plasmas 2, 3933. 

Lopreore, C. L., and R. E. Wyatt, 1999, Phys. Rev. Lett. 82, 
5190. 

Lundin, J., J. Zamanian, M. Marklund, and G. Brodin, 2007, 

Phys. Plasmas 14, 062112. 
Maafa, N., 1993, Physica Scripta 48, 351. 
Madelung, E., 1926, Z. Phys. 40, 32. 



24 



Maier, S. A., 2007, Plasmonics (Springer, New York). 
Malkin, V. M., N. J. Fisch, and J. S. Wurtele, 2007, Phys. 

Rev. E 75, 026404. 
Manfredi, G., 2005, Fields Inst. Commun. 46, 263. 
Manfredi, G., and F. Haas, 2001, Phys. Rev. B 64, 075316. 
Manfredi, G., and P.-A. Hervieux, 2007, Appl. Phys. Lett. 

91, 061108. 

Marburger, J. H., and R. F. Tooper, 1975, Phys. Rev. Lett. 
35, 1001. 

Marklund, M., and P. K. Shukla, 2006, Rev. Mod. Phys. 78, 
591. 

Marklund, M., and G. Brodin, 2007, Phys. Rev. Lett. 98, 
025001. 

Marklund, M., G. Brodin, L. Stenflo, and C. S. Liu, 2008, 

Europhys. Lett. 84, 17006. 
Markowich, P. A., et al, 1990, Semiconductor Equations 

(Springer, Berlin). 
McKinstrie, C. J., and R. Bingham, 1992, Phys. Fluids B 4, 

2626. 

Masood, W., B. Eliasson, and P. K. Shukla, Phys. Rev. E 81, 
066401 (2010). 

Mayor, F. S., A. Askar, and H. A. Rabitz, 1999, J. Chem. 

Phys. Ill, 2423. 
Melrose, D. B., 2008, Quantum Plasmadynamics: Unmagne- 

tized Plasmas (Springer, New York). Lecture Notes Phys. 

735. 

Melrose, D. B., and A. Mushtaq, 2009, Phys. Plasmas 16, 
094508. 

Mendonga, J. T., 2001, Theory of Photon Acceleration (Insti- 
tute of Physics, Bristol). 

Misra, A. P., 2007, Phys. Plasmas 14, 064501. 

Misra, A. P., D. Ghosh, and A. R. Chowdhury, 2008, Phys. 
Lett. A 372, 1469. 

Misra, A. P., 2009, Phys. Plasmas 16, 033702. 

Misra, A. P., and S. Samanta, 2010, Phys. Rev. E 82, 037401. 

Mithen, J. P., J. Daligault, and G. Gregori, 2011, Phys. Rev. 
E 83, 015401(R). 

Moyal, J. E., 1949, Proc. Cambridge Philos. Soc. 45, 99. 

Murillo, M. S., 2010, Phys. Rev. E 81, 036403 (2010). 

Murtaza, G. M., and P. K. Shukla, 1984, J. Plasma Phys. 31, 
423. 

Mushtaq, A., and D. B. Melrose, 2009, Phys. Plasmas 16, 
102110. 

National Research Council, 1995, Plasma Science: From Fun- 
damental Research to Technological Applications (National 
Academy Press, Washington DC). 

Neumayer, P., C. Fortmann, T. Dopner, et al., 2010, Phys. 
Rev. Lett. 105, 075003. 

Norreys, P. A., F. N. Beg, Y. Sentoku, et al, 2009, Phys. 
Plasmas 16, 041002. 

Oberman, C, and A. Ron, 1963, Phys. Rev. 130, 1291. 

Opher, M., et al, 2001, Phys. Plasmas 8, 2454. 

Ozbay, E., 2006, Science 311, 189. 

Pines, D., 1961, J. Nucl. Energy: Part C: Plasma Phys. 2, 5. 
Pines, D., 1983, Elementary Excitations in Solids (Benjamin, 

Massachusetts) . 
Pines, D., and P. Nozieres, 1989, The Theory of Quantum 

Liquids (Benjamin, New York). 
Redmer, R., and G. Ropke, 2010, Contrib. Plasma Phys. 50, 

970. 

Runge, E., and E. K. U. Gross, 1984, Phys. Rev. Lett. 52, 
997. 

Salamin Y. A., et al, 2006, Phys. Rep. 427, 41. 

Saleem, H., A. Ahmad, and S. A. Khan, 2008, Phys. Plasmas 



15, 014503. 
Scott R. K., 2007, Phys. Rev. E 75 046301. 
Serbeto, A., J. T. Mendonga, K. H. Tsui, and R. Bonifacio, 

2008, Phys. Plasmas 15, 013110). 

Serbeto, A., L. F. Monteiro, K. H. Tsui, and J. T. Mendonga, 

2009, Plasma Phys. Control. Fusion 51, 124024. 
Scrimaa, O. T., J. Javanainen, and S. Varro, 1986, Phys. Rev. 

A 33, 2913. 

Shaikh, D., and P. K. Shukla, 2007, Phys. Rev. Lett. 99, 
125002. 

Shaikh, D., and P. K. Shukla, 2008, New J. Phys. 10, 083007. 

Shapiro, S. L., and S. L. Teukolsky, 1983, Black Holes, White 
Dwarfs, and Neutron Stars: The Physics of Compact Ob- 
jects (John Wiley & Sons, New York). 

Sharma, R. P., and P. K. Shukla, 1983, Phys. Fluids 26, 83. 

Shukla, P. K., 2006, Phys. Lett. A 352, 242. 

Shukla, P. K., 2007, Phys. Lett. A 369, 312. 

Shukla, P. K., 2009, Nature Phys. 5, 92. 

Shukla, P. K., 2010, Phys. Lett A 374, 2897. 

Shukla, P. K., and B. Eliasson, 2006, Phys. Rev. Lett. 96, 
245001. 

Shukla, P. K., and B. Eliasson, 2007, Phys. Rev. Lett. 99, 
096401. 

Shukla, P. K., and B. Eliasson, 2008a Phys. Lett. A 372, 
2897. 

Shukla, P. K., and B. Eliasson, 2008b, Phys. Rev. Lett. 100, 
036801. 

Shukla, P. K., and B. Eliasson, 2010, Phys. Usp. 53, 51. 
Shukla, P. K., and L. Stenflo, 2006, Phys. Plasmas 13, 044505. 
Shukla, P. K., and M. Y. Yu, 1984, Plasma Phys. Controll. 
Fusion 26, 841. 

Shukla, P. K., M. Y. Yu, H. U. Rahman, and K. H. Spatschek, 
1981, Phys. Rev. A 23, 321. 

Shukla, P. K., N. N. Rao, M. Y. Yu, and N. L. Tsintsadze, 
1986, Phys. Rep. 138, 1. 

Shokri, B., and A. A. Rukhadze, 1999, Phys. Plasmas 6, 3450. 

Shpatakovskaya, G., 2006, J. Exp. Teor. Phys. 102, 466. 

Silin, V. P., and A. A. Rukhadze, 1961, Electromagnetic Prop- 
erties of Plasmas and Plasma-like Media (Gosatomizdat, 
Moscow) . 

Simpson, G., C. Sulem, and P. L. Sulem, 2009, Phys. Rev. E 
80, 056405. 

Son, S., and N. J. Fisch, 2005, Phys. Rev. Lett. 95, 225002. 
Steinberg, M., 2000, Thermodynamics and Kinetics of a Mag- 
netized Quantum Plasma (Logos, Berlin). 
Stenflo, L., and P. K. Shukla, 2009, in From Leonardo to 

ITER: Nonlinear and Coherence Aspects, Ed. J. Weiland. 

(AIP Conf. Proc, Vol 1177, New York), pp. 4-9. 
Stockman, M. I., 2011, Phys. Today 64, 39. 
Stratonovich, R. L., 1956, Sov. Phys. D 1, 414. 
Takabayasi, T., 1953, Prog. Theor. Phys. 9, 187. 
Takabayasi, T., 1955, Prog. Theor. Phys. 14, 283. 
Thiele, R., T. Bornath, C. Fortmann, A. H611, R. Redmer, 

H. Reinholz, G. Ropke, A. Wierling, S. H. Glenzer, and G. 

Gregori, 2008, Phys. Rev. E 78, 026411. 
Tsintsadze, L. N., 2010, in New Frontiers in Advanced Plasma 

Physics, Eds. B. Eliasson and P. K. Shukla, (AIP Conf. 

Proc. #1306, AIP, New York), pp. 89-102. 
Tsintsadze, N. L., and L. N. Tsintsadze, 2009, Europhys. Lett. 

88, 35001. 

Watanabe, H., 1956, J. Phys. Soc. Jpn. 11, 112. 
Wigner, E., 1932, Phys. Rev. 40, 749. 
Wilhelm, H. E., 1971, Z. Physik 241, 1. 

Wyatt, R. E., 2005, Quantum Dynamics with Trajectories: 



25 



Introduction to Quantum Hydrodynamics (Springer Sci- 
ence, New York). 
Xia, S. X., W. C. Hua, and G. Feng, 2010, Chin. Phys. Lett. 
27, 025204. 

Yu, M. Y., et al, 1974, Z. Naturforsch 29a, 1736. 



Zamanian, J., M. Marklund, and G. Brodin, 2010, New J. 

Phys. 12, 043019. 
Zhang, W., and R. Balescu, 1988, J. Plasma Phys. 40, 199. 
Zhu, J., and P. Ji, 2010, Phys. Rev. E 81, 036406. 



